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ABSTRACT 


Nickel,  cadmium,  and  chromium  are  toxic  industrial  chemicals  with  an  exposure 
risk,  found  in  military,  occupational,  and  environmental  settings.  While  the  substances 
are  known  to  have  adverse  health  effects,  the  exact  mechanisms  of  toxicity  remain 
unclear  and  a  paucity  of  biomarkers  of  exposure  and  effect  exist.  To  identify  candidate 
biomarkers  and  to  elucidate  the  mechanisms  of  toxicity  of  Ni,  Cr,  and  Cd,  H4-II-E-C3 
and  MHICI  rat  liver-derived  cell  lines  were  treated  with  various  concentrations  of  each 
metal  and  gene  expression  patterns  were  determined  through  the  use  of  an  Affymetrix 
microarray  and  analyzed  using  bioinformatic  tools  including  Ingenuity  Pathway  Analysis. 
A  total  of  992  probe  sets  were  differentially  expressed  as  a  result  of  exposure  to  nickel, 
chromium,  and/or  cadmium,  246  of  which  may  be  further  investigated  as  candidate 
biomarkers.  The  modulated  genes  were  involved  in  biological  processes  such  as  the 
oxidative  stress  response,  apoptosis,  cell  cycle  regulation,  and  hypoxia. 


INTRODUCTION 


In  the  line  of  duty  Ameriean  soldiers  are  not  only  in  danger  from  hostile  aetion, 
but  also  from  the  environment  surrounding  them.  One  sueh  risk  to  the  warfighter,  as  well 
as  to  the  general  publie,  is  the  exposure  to  Toxie  Industrial  Chemieals  (TICs)  and  Toxie 
Industrial  Materials  (TIMs).  Diagnosis  and  treatment  of  an  exposure  is  often  diffieult  as 
these  substanees  ean  have  an  irregular  distribution  in  the  environment  and  individuals 
may  respond  differently  to  them.  They  are  deemed  as  large  of  a  threat  as  chemieal 
warfare  agents  sinee  they  are  readily  available,  produeed  in  large  amounts,  and  may 
similarly  degrade  mission  performanee  (MeKone  et  al.  2000). 

Reeently,  soldiers  in  the  field  have  been  exposed  to  TICs  and  TIMs  that  were  both 
deliberately  spread  and  eneountered  in  the  environment  as  a  pollutant.  In  2003,  National 
Guard  troops  in  Iraq  were  exposed  to  hexavalent  chromium  while  protecting  contractors 
who  were  repairing  the  Qarmat  Ali  Water  Treatment  Plant  in  Basrah,  Iraq  (The 
Associated  Press  2009).  The  burn  pit  at  Joint  Base  Balad  has  also  received  considerable 
recent  media  attention  when  soldiers  reported  illnesses  they  attributed  to  inhaling  toxic 
smoke  emanating  from  the  burn  pit  (U.S.  Army  Center  for  Health  Promotion  and 
Preventive  Medicine  2008). 

The  general  public  is  also  at  risk  of  exposure  to  TICs  and  TIMs  in  the 
environment  or  as  a  result  of  occupational  hazards.  From  1980  to  1983  a  national  survey 
conducted  by  the  National  Institute  for  Occupational  Safety  and  Health  estimated  that 
727,240  workers  were  potentially  exposed  to  nickel  metal,  nickel  alloys,  or  nickel 
compounds  (ATSDR  2005a).  Additionally,  approximately  9,000  hazardous  substance 
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releases  were  reported  to  oecur  annually  in  the  15  states  taking  part  in  the  CDC 
Hazardous  Substanees  Emergeney  Events  Surveillanee  system  (ATSDR  2005b). 

The  heavy  metals  niekel  (Ni),  eadmium  (Cd),  and  ehromium  (Cr)  are  examples  of 
TICs  that  have  a  military,  oecupational,  and  environmental  relevance.  Nickel  is  used  in 
producing  batteries  and  stainless  steel,  creating  an  occupational  and  environmental 
hazard  as  discussed  above.  It  has  also  been  considered  by  the  military  for  use  as  an  alloy 
with  tungsten  and  cobalt  as  an  environmentally  friendly  replacement  for  depleted 
uranium  munitions,  but  was  dismissed  when  the  alloy  was  shown  to  be  extremely 
carcinogenic  (Kalinich  et  al.  2005).  Chromium  is  often  used  as  a  pigment  or  as  an  anti¬ 
corrosive,  making  it  a  potential  pollutant  to  the  environment  or  soldiers,  such  as  in  the 
case  at  the  Qarmat  All  water  treatment  plant  (Barceloux  1999).  Cadmium  is  frequently 
encountered  as  an  occupational  hazard  since  it  is  commonly  used  in  pigments,  batteries, 
and  welding,  a  ubiquitous  practice  in  both  the  military  and  industry  (Meo  and  Al-Khlaiwi 
2003). 

While  Ni,  Cd,  and  Cr  are  widely  used  and  have  been  well  characterized,  few 
biomarkers  currently  exist  and  their  exact  mechanisms  of  toxicity  remain  unclear.  All 
three  chemicals  are  thought  to  cause  oxidative  stress,  which  is  capable  of  inducing  DNA 
damage  and  protein  degradation  (Beyersmann  and  Hartwig  2008).  Nickel  and  chromium 
undergo  Eenton  type  reactions  forming  reactive  oxygen  species,  while  cadmium  is 
thought  to  cause  oxidative  stress  through  the  inhibition  of  antioxidant  enzymes 
(Beyersmann  and  Hartwig  2008,  Stohs  et  al.  2001).  Chromium  is  the  only  metal  of  the 
three  that  is  able  to  directly  interact  with  DNA,  forming  Cr-DNA  adducts  and  causing 
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DNA  damage,  while  niekel  and  eadmium  are  thought  to  eause  DNA  damage  through  the 
inhibition  of  repair  enzymes  (Beyersmann  and  Hartwig  2008).  Niekel  is  the  only  metal  of 
the  three  that  has  been  shown  to  mimie  hypoxia,  up-regulating  genes  involved  in  many 
eellular  proeesses  sueh  as  glyeolysis  and  apoptosis  through  the  stabilization  of  the 
transeription  faetor  HIF-la  (Salnikow  et  al.  2004).  Ni,  Cd,  and  Cr  are  also  thought  to 
eause  deregulation  of  eell  proliferation  through  various  signaling  pathways  and 
transeription  faetors,  perhaps  due  to  reaetive  oxygen  speeies  (ROS),  although  the 
aetivation  of  these  pathways  is  poorly  understood  (Beyersmann  and  Hartwig  2008). 

An  individual’s  exposure  to  these  toxie  ehemieals  is  often  diffieult  to  deteet  and 
treat  sinee  responses  and  elinieal  signs  to  these  substanees  vary  greatly;  however 
biomarkers  may  serve  as  diagnostie  tools  to  identify  an  exposure.  The  FDA  defines 
biomarkers  as  eharacteristies  that  are  objeetively  measured  and  evaluated  as  indicators  of 
normal  biologic  processes,  pathogenic  processes,  or  pharmacologic  responses  to  a 
therapeutic  intervention  (U.S. Department  of  Health  and  Human  Services,  Food  and  Drug 
Administration,  Center  for  Drug  Evaluation  and  Research,  &  Center  for  Biologies 
Evaluation  and  Research  2008).  A  chemical  or  its  metabolite  in  a  person’s  body  fluid 
may  serve  as  a  biomarker,  or  it  may  be  a  gene  expression  change  that  results  from  an 
exposure  to  a  substance.  A  biomarker  may  not  necessarily  be  a  single  characteristic,  but 
may  be  a  group  of  them  forming  a  signature  or  fingerprint  of  response.  Unfortunately, 
few  biomarkers  of  effect,  such  as  enzymes,  intermediates,  or  breakdown  products,  have 
been  identified  and  validated  for  use  in  humans  due  to  the  previous  lack  of  appropriate 
analytical  methods  and  computational  power.  Biomarker  discovery  is  a  difficult  process 
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due  to  the  lengthy  path  from  eandidate  diseovery  to  elinical  assay  and  the  laek  of  a 
rigorous  and  eomprehensive  proeess  for  biomarker  development  (Dalmasso  2008). 

The  advaneing  field  of  toxieogenomies  has  a  range  of  new  teehnologies,  sueh  as 
mieroarray  teehnology,  data  analysis,  and  bioinformatie  proeesses  and  tools,  whieh  may 
aid  in  the  diseovery  of  new  biomarkers.  Toxieogenomies  evaluates  ehanges  in  gene 
expression  in  response  to  environmental  agents  or  toxieants,  allowing  for  genome-wide 
assessment  of  gene  expression  ehanges  resulting  from  sueh  exposures.  DNA  mieroarrays 
measure  the  abundance  of  thousands  of  gene  transcripts  simultaneously,  providing  a 
snapshot  of  what  is  occurring  inside  a  cell  at  a  given  moment.  Statistical  and 
normalization  methods  have  been  developed  to  handle  differences  between  samples  and 
the  large  amounts  of  data  generated  by  a  mieroarray  experiment.  Perhaps  the  greatest 
challenge  of  mieroarray  technology  is  assigning  biological  meaning  to  the  large  amount 
of  data  generated.  Many  different  bioinformatie  programs,  applications,  and  methods 
have  been  developed  to  aid  in  deciphering  mieroarray  results,  including  programs  to 
cluster  similarly  behaving  genes  and  to  provide  insight  into  function  and  pathway 
analysis.  Therefore,  toxieogenomies  provides  tools  to  analyze  gene  expression  changes 
that  can  provide  possible  targets  for  genomic  biomarkers  as  well  as  understanding  the 
mechanisms  of  toxicity. 

In  order  to  identify  candidate  biomarkers  and  to  elucidate  the  mechanisms  of 
toxicity  of  Ni,  Cr  (VI),  and  Cd,  toxicogenomic  and  bioinformatie  methods  were  used  to 
investigate  gene  expression  patterns  of  rat  liver  based  in  vitro  systems.  H4-II-E-C3 
(H4IIE)  cells  were  selected  as  they  are  well  characterized  and  metabolically  active  liver 
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models  (Michels  et  al.  2006).  MHICI  cells,  another  well  characterized  and  liver-like  cell 
line  (Odashima  et  al.  1979),  were  then  used  to  confirm  the  results  were  reflective  of 
hepatocyte-like  function.  MHlCl  and  H4-II-E-C3  (H4IIE)  cells  were  exposed  to  nickel 
(II)  chloride  (NiCb),  cadmium  chloride  (CdCb),  and  sodium  dichromate  (Na2Cr207). 
Affymetrix  microarrays  were  used  to  evaluate  the  quantity  of  specific  mRNA  sequences 
present  in  the  samples.  Statistical  software  tools  were  used  to  determine  whether 
significant  changes  occurred  with  treatment  to  the  various  toxicants.  Bioinformatics 
software  and  applications  were  applied  to  decipher  the  biological  meaning  of  the  gene 
expression  patterns  to  aid  in  determining  the  mechanisms  of  toxicity  and  identifying 
candidate  biomarkers.  Differentially  expressed  genes  involved  in  the  known  mechanisms 
for  the  three  metals,  such  as  apoptosis,  cell  cycle  control,  DNA  damage,  and  hypoxic 
response  were  identified. 
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MATERIALS  AND  METHODS 


Cell  Culture  Conditions  and  Exposures 

MHICI  and  H4IIE  cells  (ATCC,  Manassas,  VA)  were  grown  to  near  eonflueney 
in  Dulbeeco's  Modified  Eagle's  Medium  (DMEM)  (Eonza,Walkersville,  MD)  with  10% 
fetal  bovine  serum  (EBS)  (Invitrogen,  Carlsbad,  CA)  and  10  mL  Glutamax  (Invitrogen)  at 
37  °C  and  5%  earbon  dioxide  in  T75  flasks.  Exposures  were  initiated  once  flasks  were  90 
+/-  10%  eonfluent  using  the  test  ehemicals  NiCl2-6H20,  CdCfi,  and  Na2Cr207  (Sigma- 
Aldrieh,  St.  Eouis,  MO).  Exposure  doses  were  chosen  at  the  no  observable  effect  level, 
20%  effeet,  and  50%  effeet  based  on  the  CellTiter-Eluor  Cell  Viability  and  CellTiter  96 
Aqueous  Non-Radioactive  Cell  Proliferation  Assays  (Promega,  Madison,  WI)  for  each 
cell  line  (Table  1).  Each  flask  was  washed  twice  with  serum  free  DMEM  to  remove 
residual  serum  eomponents,  as  the  media  will  be  used  in  a  parallel  proteomie  study,  with 
a  5  minute  ineubation  between  washes.  Eifteen  mL  of  serum  free  DMEM  containing  the 
proper  coneentration  of  toxieant  was  added  to  each  flask  for  24  hours.  Eour  biological 
replicates  were  performed  for  eaeh  condition,  including  a  control  containing  only  DMEM 
with  10%  EBS  and  10  mL  Glutamax.  The  medium  with  the  toxieant  was  removed  from 
eaeh  flask  and  flash  frozen  for  use  in  future  projects  after  the  24  hours  exposure  period. 

RNA  Extraction 

The  cells  were  manually  scraped  from  the  surfaee  of  the  flask  and  homogenized 
using  a  Dounee  homogenizer  (Dounce  and  Lan  1943).  Total  RNA  was  extracted  using 
Trizol  solution  per  the  manufacturer’s  instructions  for  RNA  Isolation  (Invitrogen).  An 
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Table  1:  Exposure  concentrations  (in  pM)  for  the  (a.)  H4IIE  and  (b.)  MHICI  cell  lines. 


a. 


Dose 

Nickel 

Chromium 

Cadmium 

Eow 

40 

0.275 

0.20 

Mid 

140 

1.000 

0.55 

High 

400 

10.000 

1.20 

b. 

Dose 

Nickel 

Chromium 

Cadmium 

Eow 

30 

15 

2 

Mid 

300 

180 

4 

High 

1800 

400 

8 

8 


RNeasy  Midi  Kit  cleanup  (Qiagen,  Germantown,  MD)  was  performed  as  per  the 
manufaeturer’s  instruetions  to  remove  residual  salts  and  organie  solvents.  RNA  quality 
and  quantity  was  determined  through  the  use  of  the  Agilent  Bioanalyzer  Series  II  RNA 
6000  Nano  LabChip  Kit  (Palo  Alto,  CA). 

Microarray  Preparation  and  Processing 

oDNA  and  labeled  eRNA  was  prepared  per  Affymetrix’s  GeneChip  Expression 
Analysis  Teehnieal  Manual  using  7.5  pg  total  RNA.  Twenty  mierograms  of  biotin- 
labeled  eRNA  was  sent  to  the  laboratory  of  Dr.  Maryanne  Vahey  at  the  Walter  Reed 
Army  Institute  of  Researeh  Vaeeine  Genomies  Laboratory  for  fragmentation, 
hybridization,  staining,  washing,  and  seanning  on  the  GeneChip  Rat  Genome  230  2.0 
Array  aeeording  to  the  manufaeturer’s  instruetions  (Affymetrix,  Santa  Clara,  CA). 

Data  Normalization  and  Statistical  Analysis 

Mieroarray  data  was  proeessed  for  baekground  adjustment,  normalization,  and 
summarization  using  the  Robust  Multi- Array  Averaging  method  (RMA)  (Irizarry  et  al. 
2003)  using  Partek  software,  version  6.4  Copyright  2009  (St.  Louis,  MO).  The 
mieroarray  data  was  visually  examined  for  outliers  using  Prineipal  Component  Analysis 
(PCA)  in  Partek  Pro,  a  method  of  redueing  the  multi-dimensional  data  sets  to  lower 
dimensions  that  ean  be  readily  visualized.  Samples  were  eompared  using  pairwise 
eorrelation  analysis  and  replieate  versus  replieate  dot  plots  of  all  probe  sets  in  Mierosoft 
Exeel  (Bellevue,  WA)  to  verify  inter-replieate  reprodueibility.  Replieates  were  aeeepted 
with  a  R  >0.95  and  no  gross  deviations  from  linearity.  If  a  sample  did  not  meet  these 


9 


criteria,  a  new  mieroarray  was  proeessed  from  the  total  RNA.  A  present,  absent,  or 
marginal  deteetion  eall  for  eaeh  probe  set  was  determined  by  the  Affymetrix  GCOS 
algorithm  and  only  probe  sets  with  a  present  detection  call  for  all  samples  in  at  least  one 
eondition  were  retained  for  analysis  (Areher  and  Reese  2010). 

The  probe  sets  that  met  the  present  deteetion  eriteria  underwent  a  3-way  ANOVA 
with  the  variables  ehemieal,  dose,  and  “batch”  between  eaeh  chemieal  treatment 
eondition  and  the  eontrol  using  Partek.  The  bateh  variable  was  ineluded  to  eontrol  for 
differences  observed  in  the  PCA  resulting  from  different  experimental  and  proeessing 
dates  (Figure  1).  Statistieal  signifieanee  was  determined  by  the  Benjamini  and  Hoehberg 
False  Diseovery  Rate  (FDR)  (Benjamini  and  Hoehberg  1995)  and  a  fold  ehange  filter  was 
applied.  Probe  sets  with  a  FDR  less  than  or  equal  to  0.001  and  a  1.8  or  greater  fold 
ehange  from  eontrol  in  at  least  one  treatment  eondition  were  retained  for  bioinformatie 
analysis. 

qPCR  Validation 

Quantitative  Polymerase  Chain  Reaetion  (qPCR)  was  used  to  validate  the 
mieroarray  results  from  the  toxieant  exposures  to  the  H4IIE  eell  line.  The  total  RNA  from 
the  four  biological  replicates  of  the  toxicant  exposures  used  for  mieroarray  analysis  were 
also  used  for  qPCR  validation  with  the  genes  and  primers  (Integrated  DNA  Teehnologies, 
Skokie,  IL)  shown  in  Table  2.  The  primers  were  designed  using  Primer  Express  software 
(Applied  Biosystems,  Eoster  City,  CA)  based  on  the  Referenee  Sequenee  mRNA  (Pruitt 
et  al.  2007).  eDNA  was  prepared  from  total  RNA  using  the  Advantage  RT-for-PCR  kit 
per  the  manufaeturer’s  instruetions  published  April  2006  (Cloneteeh,  Mountain  View, 
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Figure  1:  Principal  component  analysis  (PCA)  of  probe  sets  present  in  all  four  replieates 
of  at  least  one  treatment  eondition  in  H4IIE  cells.  PCA  is  a  data  reduetion  tool  used  to 
visualize  the  variability  in  an  experiment.  In  this  ease,  eaeh  data  point  represents  one 
microarray.  This  plot  illustrates  a  batch  effect  as  the  microarrays  form  two  distinct  groups 
based  on  the  dates  of  the  experiments. 
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Table  2:  Genes  and  Primers  used  for  qPCR  validation 


Gene 

Role  in  Cell 

Forward  Primer 

Reverse  Primer 

Actb 

Gapdh 

Adm 

Regulation  of  actin  eytoskeleton 

Glycolysis 

Oxidative  stress  response 

CTGTGTGGATTGGTGGCTCTATC 

GGATACTGAGAGCAAGAGAGAGGC 

CTAGCGAGGAAAAGTGCAATGC 

AGAAAGGGTGTAAAACGCAGCT 

GATGGTATTCGAGAGAAGGGAGG 

TCGCTTTAAAATGCCCTTAATACA 

Atf 

Cdt 

Fabpl 

Stress  response,  cell  cycle  control 

DNA  replication 

Oxidative  stress  response 

CACAACATTGGCGTGATTTTTT 

TTAAGCTTCCCTGTCTGCATCA 

CTATGGGTCCAAGGTGATCCAC 

CTCGGTCAGCACAGAGTAGCAC 

TCCAGAGTCTACGTCCCCTATAGC 

AATTCAGTCACGGACTTTATGCCT 

Gadd45 

DNA  repair,  cell  cycle  arrest 

GCAGAGCTGTTGCTACTGGAGAA 

TTCCCGGCAAAAACAAATAAGT 

Slc39al0 


Metal  ion  transporter 


CTACACCGGTACCATAGCTGCTT 


ACGTCTTACACGTCACTGCACC 


CA).  The  Applied  Biosystems  SYBR  Green  Master  Mix  was  used  in  a  50  pi  qPCR 
reaction  with  2  pi  of  cDNA  template  and  a  2.5  pM  final  concentration  of  each  primer  for 
a  total  of  8  genes,  2  of  which  served  as  the  endogenous  controls.  A  DNA  Opticon  2  (Bio- 
Rad,  Hercules,  CA)  was  used  for  thermal  cycling  and  fluorescence  detection  using  the 
following  scheme:  95  °C  for  10  minutes  followed  by  40  cycles  of:  95  °C  for  15  seconds, 
60  °C  for  1  minute,  and  a  fluorescence  signal  read.  Relative  fold  change  was  determined 
using  the  comparative  Ct  method  (Livak  and  Schmittgen  2001)  and  values  from  the  four 
biological  replicates  were  averaged.  The  microarray  results  were  compared  to  the  qPCR 
results  using  Pearson’s  product-moment  correlation  coefficient. 

Bioinformatics 

Multiple  exploratory  software  packages  were  used  for  bioinformatic  analysis  of 
statistically  significant  changing  genes  identified  in  the  H4IIE  cell  line.  The  probe  sets 
determined  to  be  modulated  by  toxicant  exposure  and  their  expression  values  were 
imported  into  Vxinsight  using  VxArrayImport  0.2.5  with  default  settings  (Sandia 
National  Laboratories,  Albuquerque,  NM)  for  cluster  analysis  (Kim  et  al.  2001).  Clusters 
were  manually  selected  by  their  natural  boundaries  using  the  terrain  view,  and  probe  sets 
in  each  cluster  were  exported  to  Excel  for  further  use. 

Ingenuity  Pathways  Analysis  (Ingenuity  Systems,  www.ingenuity.com,  analysis 
date  2009-11-09)  was  used  to  explore  the  biological  implications  of  the  data.  Core 
analysis  were  performed  on  the  data  using  the  Rat  Genome  230  2.0  Array  as  the  reference 
set  with  all  other  default  settings  selected.  Networks,  canonical  pathways,  and  biological 


13 


functions  that  were  statistically  significant,  as  determined  by  the  software,  were 
examined. 


Gene  Ontology  Tree  Maehine  (Zhang  et  al.  2004)  was  also  used  as  a  data 
exploration  tool.  The  probe  sets  from  each  Vxinsight  cluster  were  imported  into  the 
program  with  Rattus  norvegicus  seleeted  as  the  species  and  “rnorvegicas_affy_rat230_2” 
seleeted  as  the  gene  ID  type.  All  probe  sets  that  had  a  present  deteetion  call  in  all  four 
replicates  of  at  least  one  treatment  condition  were  uploaded  for  a  reference  set  with 
“morvegicas_affy_rat230_2”  selected  as  the  gene  ID  type.  The  hypergeometrie  method 
was  used  to  determine  statistical  significance  with  a  Benjamini  and  Hoehberg  multiple 
test  adjustment.  A  signifieance  level  of  0.05  was  used,  and  a  minimum  of  2  genes  present 
in  a  category  was  required  to  test  the  significanee  of  that  category.  The  Directed  Acyclic 
Graph  (DAG)  was  examined  for  statistically  significant  enriched  gene  ontology  (GO) 
terms. 


IPA  and  Gene  Ontology  Tree  Maehine  (GOTM)  only  identified  biological 
functions  of  a  small  subset  of  the  differentially  expressed  genes;  therefore,  a  manual 
“binning”  method  was  devised.  This  scheme  identified  the  major  biological  processes 
that  were  modulated  by  treatment  with  the  toxieants  by  assigning  probe  sets  to  groups,  or 
bins,  based  on  GO  categories.  Seven  bin  categories  were  ereated  based  on  known  effects 
of  nickel,  chromium,  and  cadmium  in  the  published  literature:  cell  cycle,  oxidative 
stress,  ion  homeostasis,  apoptosis,  energy  regulation,  hypoxic  response,  and  DNA 
damage,  replication,  and  repair.  Eaeh  bin  was  compromised  of  multiple,  related  GO  terms 
based  on  the  GO  biologieal  process  terms  provided  by  Affymetrix  in  the  annotation  file 
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for  the  Rat  Genome  230  2.0  Array  (build  29,  2009-7-13).  The  GO  terms  found  in  eaeh 
bin  ean  be  found  in  Appendix  1 .  Probe  sets  were  then  assigned  to  a  bin  if  the  GO  term 
assoeiated  with  that  probe  set  was  also  eontained  in  that  partieular  bin.  A  ehi-squared  test 
was  used  to  determine  statistieal  signifieanee  of  enriehed  bins  for  eaeh  Vxinsight  eluster. 
Expeeted  values  for  eaeh  bin  in  the  ehi-squared  test  were  ealeulated  by  the  ratio  of  probe 
sets  in  a  bin  to  the  total  number  of  probe  sets  in  a  data  set  eontaining  all  probe  sets  with 
present  deteetion  ealls  in  all  replieates  of  one  treatment  eondition  that  eontained  at  least 
one  biologieal  proeess  GO  term  in  the  Affymetrix  annotation  file.  Probe  sets  that  did  not 
eontain  any  biologieal  proeess  annotation  were  not  eonsidered  for  signifieanoe  testing. 
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RESULTS 


Data  Normalization  and  Statistical  Analysis 

The  .cel  file  of  the  microarray  data  was  RMA  preprocessed  and  outliers  were 
identified  using  a  PCA  plot,  multivariate  correlation,  and  replicate  vs.  replicate  dot  plots. 
Probe  sets  that  were  not  detected  above  the  background  level  were  removed  from  the 
analysis,  leaving  16,31 1  out  of  the  31,099  possible  probe  sets  for  further  analysis.  A  three 
way  ANOVA  (chemical,  dose,  and  batch)  was  run,  the  FDR  computed,  and  the  fold 
change  was  calculated  for  each  treatment  compared  to  the  control.  A  probe  set  with  an 
FDR  equal  to  or  less  than  0.001  and  a  fold  change  of  at  least  1.8  was  considered  to  be 
differentially  expressed.  The  analysis  identified  430  probe  sets  in  the  nickel,  456  probe 
sets  in  the  chromium,  and  288  probe  sets  in  the  cadmium  treatment  groups  as  modulated. 
Since  many  of  the  probe  sets  were  statistically  significant  in  more  than  one  chemical,  a 
total  of  992  probe  sets  overall  were  identified  as  differentially  expressed  (Figure  2  and 
Appendix  2-8).  Hierarchical  clustering  of  differentially  expressed  genes  using  Euclidean 
distance  for  row  dissimilarity  produced  a  heat  map  providing  a  visualization  of  the  probe 
sets’  expression  patterns  across  all  treatment  conditions  (Figure  3).  The  hierarchical 
clustering  showed  that  some  probe  sets  were  responding  to  only  a  single  metal  while 
others  were  responding  to  two  or  three,  and  that  response  occurred  in  a  dose  dependant 
manner. 
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Figure  2:  Venn  diagram  showing  the  distribution  of  the  992  differentially  expressed 
(FDR<0.001  and  1.8  fold  change)  probe  sets  in  the  H4IIE  cell  line  among  the  exposures 
to  nickel,  chromium,  or  cadmium. 
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Figure  3:  Hierarchical  clustering  of  992  differentially  expressed  probe  sets  showing 
some  were  responding  to  only  a  single  metal  while  others  were  responding  to  two  or  all 
three,  and  that  the  response  did  display  a  dose  dependence.  Each  column  represents  a 
treatment  condition  and  each  row  represents  the  log2  ratio  of  change  from  the  control 
sample  of  an  individual  probe  set. 
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qPCR  Validation 


Validation  of  the  microarray  results  was  conducted  by  qPCR  analysis.  Six  genes 
that  were  differentially  expressed  by  nickel,  chromium,  and/or  cadmium  were  chosen  for 
validation.  Care  was  taken  to  choose  genes  that  were  over  expressed  or  repressed  by  each 
of  the  three  chemicals.  Overall,  a  successful  validation  was  considered  to  be  at  least  a 
1.8-fold  change  in  the  same  direction  in  the  qPCR  and  microarray  experiments. 

The  comparative  Ct  method  was  used  to  determine  fold  change  differences  using 
Actb  and  Gapdh  mRNA  as  endogenous  controls.  The  qPCR  results  of  the  six  genes  used 
for  validation,  Gadd45,  Slc39al0,  Cdt,  Atf,  Fabpl,  and  Adm  successfully  validated  the 
microarray  results,  as  all  were  changing  by  at  least  two  fold  in  the  same  direction  (Figure 
4).  Pearson’s  correlation  was  also  calculated  between  the  qPCR  and  microarray  results. 
The  fold  change  of  all  genes  between  the  qPCR  and  microarray  results  were  highly 
correlated  with  r  values  for  Gadd45,  Fabpl,  Slc39al0,  Cdt,  Atf,  and  Adm  of  0.95,  0.92, 
0.96,  0.99,  0.91  and  0.97  respectively. 

Cluster  and  Pathway  Analysis  of  Microarray  Data 

The  992  probe  sets  were  imported  into  Vxinsight  to  identify  similarly  behaving 
and  presumptively  co-regulated  probe  sets  through  cluster  analysis  (Kim  et  al.  2001).  The 
cluster  analysis  groups  together  genes  with  similar  regulation  or  function  which  can  be 
explored  using  pathway  analysis  tools.  Vxinsight  uses  a  force  directed  placement 
algorithm  to  move  similar  items  closer  together  while  simultaneously  pushing  dissimilar 
objects  away  from  each  other  and  then  displays  the  relationships  on  a  3D  terrain-like  map 
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Figure  4:  Gene  expression  ehanges  detected  by  qPCR  and  microarray  are  highly 
correlated.  The  results  are  displayed  as  log2  ratio  of  change  from  control. 
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(Boyack  et  al.  2002).  Four  clusters  were  identified  and  the  probe  sets  of  each  were 
highlighted  in  either  blue,  white,  yellow,  or  green  (Figure  5  and  Appendix  2).  Eaeh  of  the 
four  clusters  represents  a  different  subset  of  probe  sets  that  were  differentially  expressed 
in  the  three  metals  (Figure  6).  The  blue  eluster  was  mainly  eomprised  of  probe  sets  that 
were  up-regulated  in  niekel  while  the  white  eluster  contained  probe  sets  that  were  mostly 
up-regulated  in  chromium  only.  The  probe  sets  in  the  green  eluster  were  up-regulated  in 
just  about  every  metal,  but  most  heavily  in  eadmium.  The  yellow  cluster  generally 
contains  probe  sets  that  are  down  regulated  in  all  three  ehemicals. 

The  probe  sets  contained  in  each  Vxinsight  cluster  were  imported  into  Ingenuity 
Pathway  Analysis  (IP A)  for  canonical  pathway  analysis  to  identify  known  pathways 
affeeted  by  treatment  to  the  toxicants  (see  Figure  7  for  a  complete  list).  A  eanonical 
pathway  in  IPA  is  a  well  eharacterized  metabolic  or  cell  signaling  pathway  that  was 
drawn  based  on  the  IPA  Knowledge  Base.  The  only  canonieal  pathway  that  was 
significantly  enriched  (p<0.05)  in  more  than  one  cluster  was  the  NRF2-mediated 
oxidative  stress  response.  The  HIFla  signaling  and  glyeolysis/glueoneogenesis  pathways 
were  also  enriehed  in  the  blue  cluster.  VDR/RXR  aetivation  and  ATM,  and  p53  signaling 
were  all  significantly  enriched  in  the  green  eluster.  In  the  white  cluster  only  the 
pyrimidine  metabolism  and  purine  metabolism  eanonieal  pathways  were  enriehed. 
Twenty-two  different  eanonieal  pathways  were  enriehed  in  the  yellow  cluster,  most  of 
which  were  involved  in  metabolism,  the  complement  system,  or  acute  phase  response. 
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Figure  5:  Cluster  analysis  of  992  differentially  expressed  probe  sets  using  Vxinsight, 
which  clusters  together  probe  sets  with  similar  expression  patterns.  Each  colored  point 
represents  one  probe  set,  and  the  height  of  each  peak  is  proportional  to  the  number  of 
data  points  beneath  it.  Probe  sets  with  similar  expression  patterns  cluster  closely  together 
while  those  that  are  different  are  further  apart.  Four  distinct  clusters  were  identified;  the 
white,  blue,  green,  and  yellow  clusters. 


22 


Blue  Cluster 


White  Cluster 


Green  Cluster 


Yellow  Cluster 


Figure  6:  Hierarchical  clustering  of  differentially  expressed  probes  sets  representing  the 
blue,  white,  green,  and  yellow  clusters  generated  by  Vxinsight  are  presented.  The  white, 
blue,  and  green  clusters  are  comprised  of  probe  sets  up-regulated  to  Cr,  Ni,  and  Cd 
respectively,  while  the  green  cluster  contains  genes  down-regulated  in  response  to  all 
three  metals.  Gene  expression  values  are  provided  as  log2  ratio  of  change. 
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Figure  7:  The  graph  displays  the  negative  log  of  Benjamini  and  Hoehberg  false  diseovery  rates  of  IPA  eanonieal  pathways  based  on 
Vxinsight  elusters  showing  pathways  enriehed  by  the  probe  sets  eontained  in  each  cluster.  The  red  line  represents  a  0.05  significance 
level  cutoff. 


The  bio  functions  of  the  four  Vxinsight  clusters  were  also  determined  using  IPA 
to  further  elucidate  mechanism  of  toxicity  (see  Figure  8  for  a  complete  list).  IPA  bio 
functions  associate  biological  functions  and  complex  interactions  from  the  Ingenuity 
Knowledge  Base  with  the  experimental  data  set.  Many  of  the  functions  involved  in  such 
processes  as  cellular  growth  and  cell  death  are  significantly  enriched  (p<0.05)  in  all 
clusters.  Five  functions,  cellular  development,  cell  cycle,  cellular  movement,  DNA 
replication,  recombination,  and  repair,  and  cellular  compromise,  are  all  enriched  in  the 
blue,  white,  and  green  cluster.  These  data  suggest  that  all  three  chemicals  are  interfering 
with  the  cell  cycle  and  causing  DNA  damage  since  the  blue,  white,  and  green  clusters 
contain  the  majority  of  up-regulated  probe  sets  across  all  chemical  treatments. 

To  determine  standardized  gene  ontology  terms  that  were  enriched,  a  process  not 
possible  using  IPA,  the  probe  sets  in  each  Vxinsight  cluster  were  imported  into  Gene 
Ontology  Tree  Machine  (GOTM).  Statistical  significance  was  determined  using  the 
hypergeometric  test  with  a  Benjamini  and  Hochberg  False  Discovery  Rate  multi  test 
correction.  GO  categories  were  deemed  to  be  significantly  enriched  at  or  below  an 
adjusted  p-value  of  0.05.  The  third  level  gene  ontology  categories  are  listed  in  Table  3. 
Using  higher  level  gene  ontology  provides  little  information,  while  a  lower  level  provides 
too  much  detail  to  list,  but  is  useful  as  an  exploratory  tool. 

IPA  and  GOTM  only  provided  function  annotation  for  a  small  amount  of  probe 
sets  identified  as  differently  expressed;  therefore,  the  probe  sets  in  each  Vxinsight  cluster 
for  which  Affymetrix  annotation  was  available  were  assigned  to  a  bin  based  on  matching 
GO  terms  (Figure  9).  The  seven  bins  were  based  on  known  effects  of  the  metals  available 
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Figure  8:  The  graph  displays  the  negative  log  of  Benjamini  and  Hoehberg  false  diseovery  rates  of  IPA  bio  function  enrichment  based 
on  Vxinsight  clusters.  The  red  line  represents  a  0.05  significance  level  cutoff. 
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Table  3:  Significant  third  level  biological  process  and  molecular  functions  of  Vxinsight 
clusters  as  determined  through  the  use  of  Gene  Ontology  Tree  Maehine  (GOTM).  The 
third  level  is  being  displayed  as  it  presents  a  usable  amount  data  without  being  too 
general. 


Cluster 

Biologieal  Proeess 

Moleeular  Eunetion 

Blue 

Cireulatory  System  Proeess 
Homeostatic  Process 
Response  to  Oxygen 

Levels 

Response  to  Nitrosative 
Stress 

E-ascorbie  Aeid  Binding 

Cation  Binding 

Identieal  Protein  Binding 

Protein  Dimerization  Aetivity 

Dioxygenase  Activity 

Oxidoreductase  Activity  Acting  on  Single  and 
Paired  Donors  with  Ineorporation  of 
Moleeular  Oxygen 

White 

Biosynthetie  Proeess 
Nitrogen  Compound 
Metabolie  Proeess 
Maeromoleeule  Metabolie 
Proeess 

Cellular  Response  to 
Stimulus 

DNA  Binding 

Purine  Nueleoside  Binding 

Hydrolase  Aetivating 

Eigase  Activity  Eorming  Nitrogen-Metal 

Bonds 

Green 

Negative  Regulation  of 
Molecular  Function 

Cell  Death 

Regulation  of  Anatomieal 
Strueture  Size 

Negative  Regulation  of 
Biologieal  Proeesses 
System  Development 

None 

Yellow 

Extraeellular  Strueture 
Organization 

Coagulation 

Digestion 

Organie  Ether  Metabolie 
Process 

Oxidation  Reduetion 
Catabolie  Proeess 
Regulation  of  Hormone 
Eevels 

Phospholipid  Binding 

Steroid  Binding 

Polysaecharide  Binding 

Sugar  Binding 

Eigase  Aetivity  forming  Carbon-Sulfur  Bonds 
Carbon-Nitrogen  Eyase  Aetivity 
Monooxygenase  Activity 

Transferase  Aetivity  Transferring  Nitrogenous 
Groups 

Peptidase  Aetivity 

Peptidase  Inhibitor  Activity 
Lipid  Transporter  Aetivity 
Heme  Transport  Aetivity 
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Figure  9:  The  graph  displays  the  bin  enrichment  significance  (-log  p-value)  by 
Vxinsight  cluster.  The  red  line  represents  a  0.05  significance  level  cutoff.  The  bin  method 
supported  the  mechanisms  of  toxicity  suggested  by  the  other  bioinformatic  applications. 
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in  the  current  literature  and  were  termed  cell  cycle,  oxidative  stress,  ion  homeostasis, 
apoptosis,  energy  regulation,  hypoxic  response,  and  DNA  damage,  replication,  and 
repair.  Statistical  significance  of  bin  enrichment  was  determined  using  a  chi  squared  test 
with  a  p-value  of  0.05  or  less  considered  significant.  In  the  blue  cluster,  the  oxidative 
stress,  ion  homeostasis,  apoptosis,  energy  regulation,  and  hypoxic  response  bins  were  all 
significantly  enriched.  In  the  white  cluster,  DNA  damage,  replication,  and  repair  and  cell 
cycle  were  the  only  bins  significantly  enriched.  Oxidative  stress,  ion  homeostasis,  and 
apoptosis  were  all  significantly  enriched  in  the  green  cluster.  The  yellow  cluster  resulted 
in  an  enrichment  of  the  DNA  damage,  replication  and  repair,  oxidative  stress,  and  cell 
cycle  bins. 

Bin  enrichment  was  also  investigated  based  on  the  individual  chemicals  to 
determine  if  any  of  the  bin  categories  were  affected  solely  by  nickel,  chromium,  or 
cadmium  exposure  (Figure  10).  A  chi-squared  test  between  a  bin  and  the  expected  value 
was  performed  to  determine  statistical  significance,  with  a  p-value  of  0.05  or  less 
considered  significant.  The  DNA  damage,  oxidative  stress,  cell  cycle,  energy  regulation, 
and  hypoxic  response  bins  were  all  enriched  in  the  probe  sets  modulated  due  to  nickel 
exposure.  The  probe  sets  with  differential  expression  due  to  chromium  treatment  resulted 
in  the  DNA  damage,  replication,  and  repair  and  oxidative  stress  bins  being  significantly 
enriched.  The  oxidative  stress,  ion  homeostasis,  apoptosis,  and  energy  regulation  bins 
were  enriched  by  the  probe  sets  modulated  by  treatment  with  cadmium.  The  bins 
significantly  enriched  are  consistent  with  the  known  effects  of  the  three  chemicals  such 
as  the  induction  ROS  causing  oxidative  stress  resulting  in  apoptosis  and  disruption  of  the 
cell  cycle. 
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Figure  10:  The  graph  shows  the  bin  enrichment  significance  by  metal  as  the  negative 
log  of  the  p-value.  The  data  supports  the  mechanisms  of  toxicity  of  the  metals  suggested 
by  the  other  bioinformatic  methods  as  well  as  those  identified  in  the  literature. 
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Candidate  Biomarker  Identification 


All  differentially  expressed  probe  sets  in  the  H4IIE  eell  line  were  initially 
eonsidered  genes  of  interest  for  further  study  as  biomarkers.  Of  the  992  differentially 
expressed  probe  sets  identified  in  the  H4IIE  eell  line,  812  were  also  differentially 
expressed  in  the  MHICI  eell  line.  Only  probe  sets  that  were  assoeiated  to  a  gene  by  the 
Affymetrix  annotation  file  were  retained,  as  this  ensures  that  only  genes  with  known 
funetion  are  being  eonsidered.  Those  genes  not  ehanging  by  2.5  fold  in  one  treatment 
eondition  were  also  discarded,  as  smaller  changes  may  be  too  small  to  serve  as  indicators 
of  exposure.  The  remaining  genes  were  then  compared  against  the  differentially 
expressed  genes  of  the  MHlCl  cell  line.  This  comparison  provided  an  opportunity  to 
validate  and  test  the  differentially  expressed  genes  identified  as  genes  of  interest  in 
another  rat  liver-derived  cell  line.  If  the  gene  was  not  changing  by  1.8  fold  in  the  MHICI 
cells  due  to  treatment  to  the  same  chemical  it  was  discarded  from  consideration  as  a 
candidate  biomarker.  Genes  that  were  not  behaving  similarly  in  both  cell  lines,  or  that 
had  conflicting  results  from  different  probe  sets  for  the  same  gene,  were  also  removed. 
This  resulted  in  9  genes  being  identified  as  candidate  biomarkers  for  all  3  chemicals,  39 
for  Cd,  99  for  Cr,  six  for  both  Cr  and  Cd,  63  for  Ni,  19  for  Ni  and  Cd,  and  1 1  for  Ni  and 
Cr  (Table  4). 
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Table  4:  Candidate  biomarker  genes  in  the  ehemicals  of  interest. 


Gene 

Chemical 

Gene 

Chemical 

Gene 

Chemical 

Aldob 

Ni,  Cd,  Cr 

ProCr 

Cd 

Ezh2 

Cr 

AtO 

Ni,  Cd,  Cr 

Ptrhl 

Cd 

E5 

Cr 

C5 

Ni,  Cd,  Cr 

RGD1560812 

Cd 

Ebxo5 

Cr 

Dde 

Ni,  Cd,  Cr 

RGD1561551 

Cd 

Eenl 

Cr 

Hmoxl 

Ni,  Cd,  Cr 

SlOOall 

Cd 

Gamt 

Cr 

Itih4 

Ni,  Cd,  Cr 

Sle3a2 

Cd 

Gea 

Cr 

Klf5 

Ni,  Cd,  Cr 

Slc41a2 

Cd 

Ghr 

Cr 

Maff 

Ni,  Cd,  Cr 

Tnfrsfl2a 

Cd 

Gins3 

Cr 

Zfand2a 

Ni,  Cd,  Cr 

TnfrsfZl 

Cd 

Hal 

Cr 

Abee2 

Cd 

Trib3 

Cd 

Hnrnpal 

Cr 

Abhd2 

Cd 

Triml6 

Cd 

Igsfll 

Cr 

AbhdS 

Cd 

Trim26 

Cd 

Inadl 

Cr 

Axudl 

Cd 

Wars 

Cd 

Iqub 

Cr 

Bag3 

Cd 

AbeaS 

Cr 

Itprl 

Cr 

Ces2 

Cd 

Adk 

Cr 

Kifl6b 

Cr 

Cnn3 

Cd 

Apob 

Cr 

Kifl8b 

Cr 

Cryab 

Cd 

Arsa 

Cr 

Klhl21 

Cr 

Ctgf 

Cd 

As3mt 

Cr 

EinVe 

Cr 

Cyp3a9 

Cd 

Bcas3 

Cr 

Erba 

Cr 

Dcbld2 

Cd 

Cacnald 

Cr 

Eyar 

Cr 

Ddit3 

Cd 

Cene2 

Cr 

Manlal 

Cr 

Dnajbl 

Cd 

Cdh2 

Cr 

Mem  10 

Cr 

Fads3 

Cd 

Cdtl 

Cr 

Mem3 

Cr 

Gadd45a 

Cd 

Chaflb 

Cr 

Mem5 

Cr 

Gadd45g 

Cd 

Cntf 

Cr 

Mem6 

Cr 

Gbp2 

Cd 

Cohhl 

Cr 

Mem7 

Cr 

Gngl2 

Cd 

Cyp4fl7 

Cr 

Mem8 

Cr 

GprlOV 

Cd 

Ddx20 

Cr 

Mpnd 

Cr 

Grhll 

Cd 

Denndla 

Cr 

Mybl2 

Cr 

Mdm2 

Cd 

Doekl 

Cr 

Mylk 

Cr 

Mlltl  1 

Cd 

Dst 

Cr 

Natl  3 

Cr 

Mobkllb 

Cd 

Elmo3 

Cr 

Ngef 

Cr 

Mxl 

Cd 

Ereebl 

Cr 

Nqol 

Cr 

Pmml 

Cd 

Exol 

Cr 

Nrbp2 

Cr 

Ppplrl5a 

Cd 

Exoe4 

Cr 

Nufipl 

Cr 
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Table  4  (Continued):  Candidate  biomarker  genes  in  the  ehemicals  of  interest. 


Gene 

Chemical 

Gene 

Chemical 

Gene 

Chemical 

Nup85 

Cr 

Traip 

Cr 

Eroll 

Ni 

Opn3 

Cr 

Uaplll 

Cr 

Errfil 

Ni 

Ore61 

Cr 

Uehl5 

Cr 

Easn 

Ni 

Osginl 

Cr 

Uox 

Cr 

Ebxo30 

Ni 

Peea 

Cr 

Uspl 

Cr 

Elotl 

Ni 

Pena 

Cr 

Weel 

Cr 

Gclc 

Ni 

Pftkl 

Cr 

Zfp367 

Cr 

Gjbl 

Ni 

PhfZOll 

Cr 

Cdr2 

Cr,  Cd 

Grbl4 

Ni 

Phldal 

Cr 

Coll8al 

Cr,  Cd 

Higdla 

Ni 

Plebl 

Cr 

Cpsl 

Cr,  Cd 

Hipl 

Ni 

Pold2 

Cr 

Lgals4 

Cr,  Cd 

Ilitm3 

Ni 

Ptges2 

Cr 

Myolg 

Cr,  Cd 

Jmjdla 

Ni 

Ptprd 

Cr 

Pig 

Cr,  Cd 

Klf6 

Ni 

Ptprg 

Cr 

Ankrd37 

Ni 

Edha 

Ni 

Qrich2 

Cr 

Apln 

Ni 

Map2kl 

Ni 

Rabgapll 

Cr 

Arlbipl 

Ni 

Mboatl 

Ni 

Rad51 

Cr 

Besll 

Ni 

Nampt 

Ni 

Rad51apl 

Cr 

BNip3 

Ni 

Nme3 

Ni 

RGD1303142 

Cr 

Btd 

Ni 

Nucb2 

Ni 

RGD 1308772 

Cr 

Cls 

Ni 

P2ry2 

Ni 

RGD 1559690 

Cr 

Cdhl7 

Ni 

P4ha2 

Ni 

Rpa2 

Cr 

Chd6 

Ni 

Pbld 

Ni 

Rrm2 

Cr 

Chsyl 

Ni 

Pdkl 

Ni 

Senp7 

Cr 

Clen4-2 

Ni 

Pgkl 

Ni 

Shank2 

Cr 

Cln3 

Ni 

Phyh 

Ni 

Skp2 

Cr 

Cptla 

Ni 

Pygl 

Ni 

Sle22al8 

Cr 

Creb311 

Ni 

Rfng 

Ni 

Sle25a21 

Cr 

Cxel9 

Ni 

RGD1561455 

Ni 

Sle48al 

Cr 

Dapll 

Ni 

RGD  1562059 

Ni 

Smoel 

Cr 

Dcxr 

Ni 

Rnf217 

Ni 

Sndl 

Cr 

Ddit4 

Ni 

Rrbpl 

Ni 

Srebfl 

Cr 

DhCr24 

Ni 

Samd4a 

Ni 

Stardl3 

Cr 

Eglnl 

Ni 

Sely 

Ni 

Tafld 

Cr 

Egln3 

Ni 

Stall 

Ni 

Tle2 

Cr 

Elovl6 

Ni 

Stbdl 

Ni 
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Table  4  (Continued):  Candidate  biomarker  genes  in  the  chemieals  of  interest. 


Gene 

Chemical 

Sultlal 

Ni 

Tes 

Ni 

Ttr 

Ni 

Txnrdl 

Ni 

Cbs 

Ni 

Ube2o 

Ni 

AmpdS 

Ni,  Cd 

Anxa2 

Ni,  Cd 

Cfi 

Ni,  Cd 

Cxcll 

Ni,  Cd 

Duspl 

Ni,  Cd 

Empl 

Ni,  Cd 

Fetub 

Ni,  Cd 

G6pe 

Ni,  Cd 

HspbS 

Ni,  Cd 

Ifrdl 

Ni,  Cd 

Jun 

Ni,  Cd 

Obfe2a 

Ni,  Cd 

Proe 

Ni,  Cd 

RGD1566118 

Ni,  Cd 

Rndl 

Ni,  Cd 

SleSOal 

Ni,  Cd 

Slc39al0 

Ni,  Cd 

Sqstml 

Ni,  Cd 

Zyx 

Ni,  Cd 

Akrlb8 

Ni,  Cr 

C9 

Ni,  Cr 

Calml4 

Ni,  Cr 

FamlVlal 

Ni,  Cr 

Fga 

Ni,  Cr 

Hmgcs2 

Ni,  Cr 

Map2k6 

Ni,  Cr 

Pah 

Ni,  Cr 

Rgn 

Ni,  Cr 

Synel 

Ni,  Cr 

Wnk4 

Ni,  Cr 
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DISCUSSION 


Nickel,  chromium,  and  cadmium  are  heavy  metals  that  have  military, 
occupational,  and  environmental  impacts.  While  the  adverse  effects  of  the  metals  have 
been  well  charaeterized,  their  exact  mechanisms  of  toxicity  have  yet  to  be  determined  and 
few  biomarkers  exist  for  the  early  detection  of  an  exposure.  The  metals  are  known  to 
cause  DNA  damage,  oxidative  stress  through  the  formation  of  reactive  oxygen  speeies, 
and/or  interfere  with  normal  cell  biological  and  physiological  processes.  Gene  expression 
studies  were  condueted,  identifying  992  Affymetrix  probe  sets  modulated  as  a  result  of 
exposure  to  niekel,  chromium,  or  cadmium.  The  use  of  bioinformatie  tools  enabled  the 
identification  of  pathways  and  functions  with  an  over-representation  of  differentially 
expressed  genes. 

qPCR  Validation 

Validation  of  microarray  data  is  a  common  practice  since  the  quality  of  the  data 
can  vary  with  the  platform  and  procedures  used.  Quantitative  real-time  PCR  (qPCR)  is  a 
technique  commonly  used  to  validate  mieroarray  results.  Correlations  between  the 
microarray  and  qPCR  data  were  ealeulated  using  Pearson’s  eorrelation  with  results 
between  0.91  and  0.97.  All  genes  used  for  validation  that  were  changing  by  at  least  two 
fold  due  to  a  ehemieal  exposure  in  the  mieroarray  results,  were  also  ehanging  by  two  fold 
in  the  same  direetion  in  the  qPCR  results.  Since  the  fold  changes  were  similar  and  the 
correlations  were  good,  this  will  be  considered  a  suceessful  validation.  The  small 
diserepancies  between  the  mieroarray  and  qPCR  results  may  be  due  to  the  different 
effieiencies  of  the  reverse  transcription  reactions,  the  differences  in  methods  including 
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the  hybridization  step  required  only  in  mieroarray  proeessing,  or  the  different 
normalization  teehniques  used  by  the  different  methods. 

Candidate  Biomarker  Identification 

Identifying  and  validating  biomarkers  for  elinieal  use  is  a  eomplex  proeess.  The 
first  step,  identifying  differentially  expressed  genes  for  further  investigation  as  eandidate 
biomarkers,  was  undertaken  as  part  of  this  projeet.  Of  the  992  probe  sets  differentially 
expressed  in  the  H4IIE  eell  line  as  a  result  of  exposure  to  at  least  one  eoneentration  of 
nickel,  cadmium,  and/or  chromium,  246  were  identified  for  further  investigation  for  use 
as  potential  biomarkers.  A  biomarker  does  not  need  to  be  a  single  piece  of  data,  but  could 
be  a  panel  that  when  taken  as  a  whole  provides  useful  information.  A  biomarker 
comprised  of  multiple  pieces  of  information  may  more  accurately  examine  biological 
processes  on  a  system-wide  scale  (Aardema  and  MacGregor  2002).  Therefore,  these 
genes  alone  may  not  serve  as  a  viable  biomarker,  but  a  group  of  them  may  serve  as  a 
fingerprint  or  signature  of  exposure.  The  differentially  expressed  genes  were  identified  in 
one  cell  line  and  then  validated  in  another  to  ensure  the  effects  being  seen  were  typical  of 
the  liver  and  likely  to  occur  in  whole  animal  systems  as  well.  Further  investigation  and 
validation  of  the  modulated  probe  sets  is  warranted  as  they  were  modulated  in  a  rodent 
cell  line  and  may  not  behave  similarly  in  vivo  in  rodents  or  humans.  To  validate  these 
gene  transcripts  for  use  as  a  biomarker  of  exposure  or  effect,  the  three  metals  would  need 
to  be  tested  in  a  rat  whole  animal  system.  The  transcripts  behaving  similarly  in  vivo 
would  also  need  to  be  tested  in  a  human  liver  cell  line  to  see  if  the  genes  are  expressed  as 
part  of  a  highly  conserved  biological  process  before  being  validated  in  humans.  Also,  it 
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would  need  to  be  determined  if  the  expression  of  the  eandidate  biomarkers  is  detectable 
by  minimally  invasive  procedures,  such  as  through  a  blood  sample,  necessitating  that  the 
protein  be  expressed  in  a  detectable  manner,  such  as  being  secreted  from  the  cell.  No 
formal  procedure  for  biomarker  discovery  exists,  but  these  differentially  expressed  genes 
serve  as  a  starting  point  and  warrant  further  investigation  and  validation. 

Chromium  Toxicity  Pathways 

Vxinsight,  a  data  mining  tool,  was  used  to  cluster  probe  sets  with  similar 
expression  patterns  among  all  the  treatment  groups.  The  white  cluster  was  comprised  of  a 
total  of  129  probe  sets,  124  of  which  were  differentially  expressed  due  to  chromium 
exposure  and  8  were  differentially  expressed  in  nickel,  5  of  which  were  unique  only  to 
nickel.  Of  the  5  probe  sets  not  modulated  in  chromium,  only  3  had  annotation  through 
Affymetrix  {Spc24,  Tcfl9,  and  Tkl).  Fifty-one  of  the  129  probes  in  the  white  cluster  had 
no  functional  annotation.  Therefore,  it  is  possible  to  assign  the  effects  observed  in  this 
cluster  as  a  result  of  exposure  to  chromium,  and  allowing  it  to  be  used  to  investigate  the 
mechanisms  of  toxicity  of  chromium. 

The  white  cluster  is  dominated  by  genes  involved  in  the  DNA  damage  response 
suggesting  a  mechanism  of  toxicity  for  chromium  (Figure  1 1).  Of  the  79  probe  sets  with 
functional  annotation  in  the  Affymetrix  annotation  file,  40  were  placed  in  the  DNA 
damage,  replication,  and  repair  bin  base  on  GO  terms.  This  is  the  most  significantly 
enriched  GO  category  of  the  entire  data  set.  Chromium  has  been  shown  to  cause  DNA 
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Figure  11:  The  presented  image  is  a  network  of  genes  from  the  white  eluster  involved  in 
the  DNA  damage  response  produeed  in  IPA.  The  network  displays  the  interaetion 
between  genes  involved  in  DNA  replication  and  repair  among  themselves  and  with  genes 
involved  in  other  biological  functions  such  as  cell  cycle  control.  Genes  with  similar 
functions  such  as  DNA  replication  initiation  and  double  strand  break  repair  are  labeled  as 
such.  The  intensity  of  the  pink  color  represents  the  magnitude  of  the  fold  change,  with  a 
deeper  the  color  representing  a  greater  magnitude  of  change.  Interactions  are  represented 
as  lines  with  the  direction  of  the  arrow  signifying  the  directionality  of  the  interaction. 
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damage  through  the  generation  of  ROS  and  by  direetly  forming  Cr-DNA  adduets 
(Beyersmann  and  Hartwig  2008). 

The  up-regulation  of  genes  involved  in  DNA  replieation  and  repair,  sueh  as  DNA 
polymerases  and  those  that  make  up  the  DNA  synthesome,  also  suggest  that  ehromium 
eauses  DNA  damage  (Bridgewater  et  al.  1998,  Dai  et  al.  2009).  The  DNA  synthesome  is 
a  multiprotein  eomplex  involved  in  DNA  replieation  eontaining  DNA  polymerases  a,  6 
{Fold),  and  s  {Pole),  DNA  primase,  topoisomerases,  proliferation  eell  nuelear  antigen 
{Pena),  replieation  faetor  C  {Rfc),  replieation  protein  A  {Rpa),  3 ’>5’  and  5 ’>3’  DNA 
helieases,  DNA  methyltransferase,  poly  (ADP)-ribose  polymerase  and  DNA  ligase  (Coll 
et  al.  1996,  Jiang  et  al.  2002).  Many  of  the  genes  eoding  for  the  proteins  in  the  DNA 
synthesome  are  differentially  expressed  in  the  white  eluster  ineluding  Pena,  Rpa,  the 
minichromosome  maintenance  complex  components  which  serve  as  helieases,  Pole,  and 
Poldl,  providing  a  method  through  which  the  cells  may  be  responding  to  the  DNA 
damage  caused  by  exposure  to  chromium. 

The  differentially  expressed  genes  that  make  up  part  of  the  synthesome  often  play 
specific  roles  in  DNA  synthesis,  suggesting  more  detailed  mechanisms  of  toxicity.  Since 
DNA  polymerases  6  and  s  are  associated  with  proof-reading  and  repair  activity,  their  up- 
regulation  may  be  in  response  to  chromium  specific  DNA  damage  (Tran  et  al.  1999). 
DNA  ligase  is  involved  in  repairing  double  strand  breaks  which  are  known  to  accumulate 
due  to  chromium  toxicity.  Some  of  the  genes  differentially  expressed  in  the  DNA 
synthesome,  such  as  the  minichromosome  maintenance  complex  (MCM)  and  origin 
recognition  complex  (ORC),  are  involved  specifically  in  the  initiation  of  DNA  synthesis 
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(Kawakami  and  Katayama  2010).  It  has  been  shown  that  Chromium-DNA  adducts  reduce 
the  ability  of  synthesome  to  initiate  replication  (Dai  et  al.  2009);  however  these  genes 
may  be  up-regulated  as  the  cell  attempts  to  repair  damage  caused  by  the  chromium. 
Additionally,  genes  involved  in  excision  repair  to  remove  DNA  adducts,  such  as  Rpa2, 
Fold,  and  Pena,  were  up-regulated  in  the  white  cluster  suggesting  that  the  direct 
interaction  of  Cr  with  DNA  and  the  formation  of  adducts  are  involved  in  Cr  toxicity. 

The  other  biological  process  in  the  binning  method  that  was  significantly  enriched 
in  the  white  cluster  was  cell  cycle.  DNA  replication  and  cell  cycle  control  are  thought  to 
be  closely  associated  (Sugimoto  et  al.  2004).  In  the  Gi  phase,  CDTl  is  recruited  to 
chromatin  bound  ORC  functioning  as  a  loader  for  the  MCM  complex,  permitting  the  cell 
cycle  to  continue  (Ishimi  1997,  Sugimoto  et  al.  2004).  Cyclin  dependent  kinases  then 
repress  the  MCM’s  helicase  activity  in  the  S  phase  (Itzhaki  et  al.  1997).  As  these  genes, 
including  Cdtl,  Ccnel  which  is  a  cyclin,  and  Skp2,  a  kinase,  are  up-regulated,  the  cells 
may  have  overcome  cell  cycle  arrest  as  a  survival  technique  for  resistance  to  the  toxic 
insult.  Therefore,  chromium  toxicity  may  interfere  with  cell  cycle  regulation  and  induce 
the  DNA  replication  machinery. 

Nickel  Toxicity  Pathways 

The  blue  cluster,  as  determined  through  Vxinsight,  contained  216  differentially 
expressed  probe  sets,  1 1 1  of  which  had  no  functional  annotation  provided  by  Affymetrix. 
Of  the  216  differentially  expressed  probe  sets  only  4  were  modulated  only  by  Cr,  190 
solely  by  Ni,  6  by  both  Ni  and  Cd,  14  by  both  Ni  and  Cr,  and  2  by  all  three  metals  (Figure 
12).  In  the  binning  methods,  the  oxidative  stress,  ion  homeostasis,  apoptosis,  energy 
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Figure  12:  Venn  diagram  of  differentially  expressed  probe  sets  in  the  blue  eluster,  green 
eluster,  and  yellow  cluster  as  determined  through  cluster  analysis  using  Vxinsight.  The 
blue  cluster  contains  mostly  probe  sets  modulated  to  exposure  to  nickel,  the  green  cluster 
contains  probe  sets  modulated  in  response  to  exposure  to  cadmium,  and  the  yellow 
cluster  is  not  specific  to  one  chemical. 
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regulation,  and  hypoxic  response  categories  were  all  enriched  by  the  probe  sets  present  in 
the  blue  cluster.  Since  all  but  four  genes  were  differentially  expressed  in  Ni,  it  is  fair  to 
consider  effects  observed  in  this  cluster  being  a  result  of  exposure  to  Ni  and  to  use  the 
blue  cluster  to  investigate  the  mechanisms  of  toxicity  of  Ni. 

Nickel  is  thought  to  cause  toxicity  by  inducing  the  formation  of  reactive  oxygen 
species  (ROS)  in  a  Fenton  type  reaction,  catalyzing  the  generation  of  hydroxyl  radicals 
from  hydrogen  peroxide  causing  damage  to  proteins  and  DNA  (Beyersmann  and  Hartwig 
2008).  The  ROS  may  also  regulate  transcription  factors  and  signal  transduction  pathways 
involved  in  apoptosis  and  cell  proliferation  (Valko  et  al.  2005).  Many  genes  involved  in 
the  oxidative  stress  response  were  up-regulated  in  the  blue  cluster,  including,  Gpxl,  Gdc, 
Eroll,  AldhSal,  Loxl2,  and  Txnrdl,  suggesting  that  ROS  play  a  role  in  nickel  toxicity 
(Figure  13).  Glutathione  peroxidase  1  {Gpxl)  is  one  of  the  most  important  antioxidant 
enzymes  involved  in  the  detoxification  of  hydrogen  peroxide  (Brigelius-Flohe  2006). 
Glutathione  is  an  antioxidant  known  to  prevent  the  nickel  induced  generation  of  ROS  and 
the  gene  responsible  for  the  first  rate-limiting  step  in  glutathione  production,  Gdc,  was 
up-regulated  in  the  blue  bin  (Das  et  al.  2008,  Mascart  and  Gottignies  1979).  Eroll  is 
involved  in  redox  homeostasis  and  is  induced  in  the  course  of  the  unfolded  protein 
response  (Pagani  et  al.  2000),  and  may  therefore,  be  up-regulated  in  response  to  damaged 
proteins  caused  by  oxidative  stress.  Aldehyde  dehydrogenase  3A1  {Aldh3al)  plays 
multiple  roles  in  protecting  against  oxidative  stress  such  as  scavenging  ROS  and  perhaps 
reducing  cell  growth  (Estey  et  al.  2007). 
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Figure  13:  Network  of  genes  from  the  blue  cluster  involved  in  the  HIF-la  and  oxidative 
stress  response  produced  in  IPA.  Twenty-seven  genes  in  this  cluster  are  known  to  be 
regulated  by  HIF-la.  The  intensity  of  the  pink  color  represents  the  magnitude  of  the  fold 
change,  with  a  deeper  the  color  representing  a  greater  magnitude  of  the  fold  change. 
Interactions  are  represented  as  lines  with  the  direction  of  the  arrow  signifying  the 
directionality  of  the  interaction. 


43 


Nickel  has  also  been  shown  to  activate  the  hypoxia-indueible  faetor-la  (HIF-la) 
which  induces  the  transeription  of  genes  involved  in  glyeolysis,  glueose  transport, 
apoptosis,  and  other  eellular  proeess  (Beyersmann  and  Hartwig  2008).  Niekel  prevents 
the  degradation  of  HIF-la  either  by  the  depletion  of  aseorbate  or  by  replaeing  iron  in  the 
hydroxylases  responsible  for  HIF-la  degradation  (Maxwell  and  Salnikow  2004, 

Salnikow  et  al.  2004).  In  the  blue  bin,  27  out  of  the  105  differentially  expressed  probe 
sets  with  annotation  were  HIF-la  regulated  according  to  IP  A,  resulting  in  a  signifieant 
enrichment  of  the  HIF-la  eanonieal  pathway  in  IP  A,  as  well  as  the  hypoxie  response  bin 
(Figure  13).  Two  funetions  regulated  by  HIF-la,  glycolysis  and  apoptosis,  were  also 
enriehed  by  the  differentially  expressed  probe  sets  in  the  blue  eluster  using  the  binning 
method.  Genes  involved  in  glyeolysis  are  up-regulated  as  the  eell  begins  to  use  glucose  as 
the  main  energy  souree  in  the  absenee  of  oxygen  or  HIF-la  aetivation  (Marin-Hernandez 
et  al.  2009,  Semenza  2007). 

As  a  result  of  mimicking  hypoxia  and  the  aetivation  of  HIF-la,  many  of  the  genes 
involved  in  aerobic  energy  produetion  would  be  expeeted  to  be  down-regulated  while 
inereasing  glyeolysis.  The  biologieal  functions  lipid  and  carbohydrate  metabolism  were 
down  regulated  overall  in  the  yellow  eluster,  a  trend  eontinued  when  only  the  niekel 
modulated  genes  were  observed.  For  example,  genes  involved  in  lipid  metabolism  such 
as  Mgll  and  Cptla  are  down  regulated  in  the  yellow  eluster  by  niekel  only. 

The  biologieal  proeesses  modulated  by  the  genes  in  the  blue  eluster  suggest  two 
main  methods  through  whieh  niekel  may  exert  its  toxieity.  Genes  involved  in  apoptosis 
and  glyeolysis,  proeesses  regulated  by  both  ROS  and  HIF-la,  are  differentially  expressed 
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in  this  cluster.  Both  these  pathways  may  contribute  to  niekel  toxieity,  but  further 
investigation  is  warranted  to  determine  the  involvement  of  eaeh. 

Cadmium  Toxicity  Pathways 

The  green  cluster,  as  determined  through  Vxinsight,  contained  a  total  of  240 
probe  sets,  only  154  of  whieh  had  functional  annotation  provided  by  Affymetrix.  Of  the 
240  probe  sets  modulated  in  this  cluster,  18  of  the  probe  sets  were  differentially 
expressed  by  all  three  chemicals,  161  by  Cd  only,  9  by  Cd  and  Cr,  24  by  Ni  and  Cd,  17 
by  Cr  only,  2  by  Ni  and  Cr,  and  9  by  only  Ni  (Figure  12).  Therefore,  taking  into  aecount 
that  only  a  few  of  these  probe  set  are  not  differentially  expressed  as  a  result  of  exposure 
to  Cd,  the  green  cluster  eould  be  seen  as  probe  sets  that  respond  to  Cd  and  be  used  to 
investigation  the  meehanisms  of  toxicity  of  Cd. 

Many  of  the  genes  in  the  green  eluster  and  differentially  expressed  due  to 
exposure  to  Cd  may  have  been  up-regulated  through  p53  aetivation  (Figure  14). 

Cadmium  has  been  shown  to  induce  apoptosis  through  the  aetivation  of  p5  3 -mediated 
signaling  (Son  et  al.  2010).  Many  genes  regulated  by  p53  and  involved  in  apoptosis,  such 
as  Casp4,  Abcbl,  and  Mdm2,  are  up-regulated  in  the  green  cluster  and  due  to  eadmium 
intoxication.  Additionally,  biological  processes  such  as  cell  eycle,  cell  death,  and  p53 
signaling  were  enriehed  by  probe  sets  in  the  green  cluster  as  observed  in  IPA  analysis  and 
the  binning  method. 

Other  genes  in  the  green  cluster  and  modulated  by  Cd,  sueh  as  heat  shock  proteins 
(HSP)  may  be  involved  in  general  stress.  Heat  shock  proteins  are  cellular  ehaperone 
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Figure  14:  Network  of  genes,  produced  with  IP  A,  differentially  expressed  in  the  green 
cluster  and  as  a  result  of  exposure  to  cadmium.  Most  of  these  genes  are  involved  in 
apoptosis  and  are  regulated  by  p53,  suggesting  that  Cd  induces  p53-mediated  apoptosis. 
The  intensity  of  the  pink  color  represents  the  magnitude  of  the  fold  change,  with  a  deeper 
the  color  representing  a  greater  magnitude  of  the  change.  Interactions  are  represented  as 
lines  with  the  direction  of  the  arrow  signifying  the  directionality  of  the  interaction. 
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proteins  involved  in  the  folding  and  degradation  of  proteins  as  well  as  apoptosis  (Bertin 
and  Averbeek  2006).  Cadmium  has  been  shown  to  up-regulate  HSP  (Eseobar  et  al.  2009) 
whieh  may  be  responding  to  the  protein  damage  caused  by  oxidative  stress.  In  the  green 
cluster,  heat  shock  proteins  such  as  Hspala,  Hspalb,  and  HspbS  were  differentially 
expressed,  suggesting  that  Cd  induced  oxidative  stress  may  cause  protein  degradation 
Effects  Common  to  Nickel,  Chromium,  and  Cadmium 

Eighteen  of  the  probe  sets  in  the  green  cluster  are  differentially  expressed  in  all 
three  chemical  exposures.  The  genes  in  this  group  are  comprised  mostly  of  those 
involved  in  general  and  oxidative  stress  or  metal  response  such  as  heme  oxygenase  1 
(Hmoxl),  metallothionein  la  and  2a  {Mtla  and Mt2a),  and  activating  transcription  factor 
3  {Atf2).  Hmoxl  is  a  ubiquitous  stress  response  protein  involved  in  reducing  the  effects  of 
oxidative  stress  and  apoptosis  (Ryter  and  Choi  2009).  Metallothioneins  are  involved  in 
the  transport  and  detoxification  of  heavy  metals  and  protection  against  oxidative  stress 
(Krizkova  et  al.  2009).  Atf3  is  a  transcription  factor  whose  expression  has  been  shown  to 
increase  in  response  to  stress  (Hai  et  al.  1999). 

Additionally,  the  yellow  cluster  contains  1 1  probe  sets  differentially  expressed 
due  to  exposure  to  all  three  of  the  chemicals  tested.  Eour  hundred-seven  probe  sets  were 
down  regulated  throughout  all  three  metals;  156  in  nickel,  265  in  chromium,  and  68  in 
cadmium,  with  some  overlap  between  the  chemicals  (Eigure  12).  Many  of  the  genes  in 
this  cluster  may  be  down  regulated  in  response  to  oxidative  stress,  as  that  bin  was 
significantly  enriched.  Eor  example,  Glutaredoxin-1  (GERX),  a  protein  whose  activity 
has  been  shown  to  be  decreased  by  both  hydrogen  peroxide  and  cadmium,  is  involved  in 
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apoptosis  and  the  oxidative  stress  response.  Also,  many  of  the  genes  in  the  yellow  eluster, 
sueh  as  Ctnnbl,  Itgbl,  eadherins  Cdh2  and  Cdhl7,  are  involved  in  eell  adhesion 
(Prozialeek  et  al.  2002).  These  genes  may  be  down  regulated  as  a  result  of  the  eells 
detaehing  from  the  surfaee  and  eaeh  other  as  a  result  of  apoptosis  and  interruption  of  the 
eell  eyele. 

Interestingly,  using  the  binning  method  for  the  yellow  eluster  resulted  in  two  bins 
being  signifieant  as  a  result  of  eontaining  fewer  genes  than  expeeted  than  by  ehanee 
alone.  This  may  be  less  of  a  result  of  biologieal  funetion  and  more  of  a  result  of  the 
method  used.  The  two  bins  that  eontained  signifieantly  fewer  genes  than  expeeted  were 
the  DNA  damage,  replieation,  and  repair  and  the  eell  eyele  bin.  It  is  possible  that  sinee 
these  two  bins,  espeeially  DNA  damage,  replieation  and  repair,  were  so  strongly  enriehed 
in  other  elusters  that  not  as  many  probe  sets  would  have  been  available  to  enrieh  those 
bins  in  the  yellow  eluster. 

Future  Directions 

While  mieroarray  analysis  suggests  toxieity  pathways  involved  in  niekel, 
ehromium,  and  eadmium  intoxieation,  hard  eonelusions  about  the  exaet  meehanisms 
eannot  be  made  without  further  investigation.  This  projeet  has  suggested  that  the  DNA 
repair  and  replieation  proeess  has  been  affeeted  by  ehromium,  perhaps  through  disruption 
of  the  DNA  synthesome  or  direetly  damaging  the  DNA.  Sinee  Cr  produees  ROS  and 
forms  Cr-DNA  adduets,  mieroarray  analysis  alone  eannot  determine  the  extent  to  whieh 
eaeh  plays  a  role  in  Cr  toxieity  and  DNA  damage.  Studies  using  a  ROS  seavenger  to 
remove  oxidative  stress  may  be  useful  in  determining  the  role  the  two  meehanisms  play 
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in  chromium  toxicity.  Similarly,  many  of  the  differentially  expressed  genes  as  a  result  of 
exposure  to  nickel  may  be  controlled  by  the  transcription  factor  HIF-la.  However, 
mieroarray  teehnology  cannot  detect  HIF-la  activation  since  HIF-la  is  regulated  at  the 
protein  level,  and  niekel  blocks  its  degradation.  A  western  blot  eould  be  used  to  measure 
the  amount  of  the  HIF-la  in  cells  exposed  to  nickel  to  see  if  the  protein  is  stabilized. 

The  extensive  data  gathered  not  only  provides  many  avenues  of  exploration,  but 
also  complieates  the  analysis  itself.  Hundreds  of  bioinformatic  tools  are  available  to 
deeipher  the  biological  meaning  of  microarray  experiments,  however  none  are  perfect.  In 
this  project  IP  A,  GOTM,  Vxinsight,  and  others  initially  used  for  initial  data  exploration 
were  used  for  analysis.  These  packages  cover  free,  web-based  applications  to 
commercial,  manually  curated  software  packages.  While  these  tools  do  aid  in  annotation 
and  exploration,  especially  when  applying  statistieal  signifieance,  nothing  replaces 
intensive  manual  annotation  and  investigation. 

Additionally,  992  probe  sets  were  identified  as  differentially  expressed;  however, 
only  roughly  a  third  of  them  had  any  functionally  annotation  through  current 
bioinformatie  methods.  The  genes  without  annotation  provide  another  opportunity  for 
further  investigation  as  they  may  play  significant,  and  currently  unknown,  roles  in  the 
cell’s  response  to  Ni,  Cr,  and  Cd.  Since  many  of  these  unknown  gene  transcripts 
clustered  with  genes  of  known  function,  they  may  have  similar  roles  or  be  regulated  by 
the  same  mechanisms.  Basic  Local  Alignment  Search  Tool  (BLAST)  could  be  used  to 
determine  if  the  sequences  of  these  probe  sets  are  homologous  to  any  known  genes  or  if 
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their  promoter  sequenees  eontain  binding  sites  for  the  same  transeription  faetors, 
suggesting  similar  eontrol  meehanisms  (Altsehul  et  al.  1990). 

Conclusion 

Niekel,  ehromium,  and  eadmium  are  heavy  metals  found  in  the  military,  industry, 
and  the  environment  that  have  adverse  health  effeets  on  individuals  exposed  to  them. 
While  the  meehanisms  of  toxieity  of  the  metals  are  not  eompletely  understood,  it  is 
thought  that  they  have  similar  effeets,  but  with  different  pathways.  In  order  to  investigate 
the  effects  of  these  metals  and  to  investigate  possible  biomarkers,  cells  were  treated  with 
nickel,  chromium,  or  cadmium  and  the  gene  expression  patterns  were  measured  using 
microarray  technology.  A  total  of  992  unique  probe  sets  were  differentially  expressed  as 
a  result  of  treatment  to  nickel,  chromium,  or  cadmium,  suggesting  unique  molecular 
mechanisms  of  action  for  each  metal. 

Genes  whose  transcripts  were  differentially  expressed  in  both  the  H4IIE  and 
MHICI  cell  line  were  identified  as  genes  of  interest  for  further  investigation  as 
biomarkers.  A  total  of  246  genes  with  known  functions  were  modulated  similarly  in  both 
cell  lines.  Identifying  the  genes  of  interest  in  one  rat  liver-derived  cell  line,  and  validating 
the  changes  in  another  rat  liver-derived  cell  line  ensured  that  the  results  were  reflective  of 
hepatocyte-like  function. 

Bioinformatic  tools  were  used  to  determine  biological  functions  affected  by 
exposure  to  nickel,  chromium,  and/or  cadmium.  The  992  differentially  expressed  probe 
sets  were  shown  to  be  involved  in  such  biological  processes  as  oxidative  stress  response. 
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interference  with  the  cell  cycle,  and  apoptosis.  All  three  metals  increased  the  expression 
of  genes  in  response  to  general  and  metal  stress,  such  as  Atfl  and  metallothioneins,  and 
oxidative  stress,  such  as  Hmox.  Chromium  is  the  only  metal  of  the  three  to  directly 
interact  and  cause  damage  to  DNA,  as  observed  with  the  heavily  enriched  DNA  damage, 
replication,  and  repair  bin  of  the  white  cluster.  Nickel  was  shown  to  induce  oxidative 
stress,  and  interfere  with  energy  production  and  induce  apoptosis  perhaps  through  HIF-la 
signaling.  Cadmium  increased  expression  of  genes  as  a  result  of  oxidative  stress  such  as 
the  HSPs  and  caused  activation  of  the  p53  pathway  inducing  cell  cycle  arrest  and 
apoptosis. 

In  conclusion,  the  gene  expression  of  two  rat  liver  hepatoma  cell  lines  were 
investigated  to  discover  candidate  biomarkers  of  exposure  and  effect  and  identify 
possible  mechanism  of  toxicity  of  nickel,  chromium,  and  cadmium.  A  total  of  992  unique 
differentially  expressed  probe  sets  were  identified,  246  of  which  are  reproducible  and 
may  be  further  investigated  as  possible  biomarkers  of  effect  for  the  chemicals.  Potential 
mechanisms  of  toxicity  with  unique  molecular  modes  of  action  were  identified  for  the 
metals  including  oxidative  stress,  DNA  damage,  and/or  hypoxia,  providing  avenues  for 
further  investigation. 
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APPENDIX  1:  GO  TERMS  FOR  BIN  CATEGORIES 


DNA  Damage,  Replication,  and  Repair 

base-excision  repair,  DNA  ligation 
cellular  response  to  DNA  damage  stimulus 
DNA  catabolic  process 
DNA  catabolic  process,  exonucleolytic 
DNA  damage  checkpoint 

DNA  damage  response,  detection  of  DNA  damage 

DNA  damage  response,  signal  transduction 

DNA  damage  response,  signal  transduction  by  p53  class  mediator 

DNA  damage  response,  signal  transduction  by  p53  class  mediator  resulting  in  cell  cycle 

arrest 

DNA  damage  response,  signal  transduction  by  p53  class  mediator  resulting  in  induction 
of  apoptosis 

DNA  damage  response,  signal  transduction  by  p53  class  mediator  resulting  in 
transcription  of  p21  class  mediator 

DNA  damage  response,  signal  transduction  resulting  in  induction  of  apoptosis 
DNA  dealkylation 

DNA  double-strand  break  processing 

DNA  duplex  unwinding 

DNA  integration 

DNA  ligation 

DNA  metabolic  process 

DNA  methylation 

DNA  methylation  on  cytosine  within  a  CG  sequence 

DNA  recombination 

DNA  repair 

DNA  replication 

DNA  replication  checkpoint 

DNA  replication  initiation 

DNA  replication  proofreading 

DNA  replication,  synthesis  of  RNA  primer 

DNA  synthesis  during  DNA  repair 

DNA  topological  change 

DNA  unwinding  during  replication 

DNA-dependent  DNA  replication 

error-prone  postreplication  DNA  repair 

intra-S  DNA  damage  checkpoint 

maintenance  of  DNA  methylation 

maintenance  of  DNA  repeat  elements 

mitochondrial  DNA  replication 

negative  regulation  of  DNA  binding 

negative  regulation  of  DNA  recombination 

negative  regulation  of  DNA  repair 
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APPENDIX  1  (CONTINUED) 


DNA  Damage,  Replication,  and  Repair  (Continued) 

negative  regulation  of  DNA  replication 

negative  regulation  of  transcription  from  RNA  polymerase  II  promoter  in  response  to 
UV-induced  DNA  damage 

nucleotide-excision  repair,  DNA  damage  removal 
nucleotide-excision  repair,  DNA  duplex  unwinding 
nucleotide-excision  repair,  DNA  gap  fdling 
nucleotide-excision  repair,  DNA  incision,  3 '-to  lesion 
nucleotide-excision  repair,  DNA  incision,  5 '-to  lesion 
positive  regulation  of  DNA  binding 

positive  regulation  of  DNA  damage  response,  signal  transduction  by  p53  class  mediator 

positive  regulation  of  DNA  ligation 

positive  regulation  of  DNA  repair 

positive  regulation  of  DNA  replication 

regulation  of  DNA  binding 

regulation  of  DNA  recombination 

regulation  of  DNA  replication 

regulation  of  DNA  replication  initiation 

response  to  DNA  damage  stimulus 

M  transition  DNA  damage  checkpoint 

M  transition  DNA  damage  checkpoint 

centromeric  heterochromatin  formation 

ATP-dependent  chromatin  remodeling 

chromatin  assembly  or  disassembly 

chromatin  modification 

chromatin  remodeling 

chromatin  silencing 

chromatin  silencing  at  telomere 

chromosome  condensation 

chromosome  organization 

chromosome  segregation 

establishment  or  maintenance  of  chromatin  architecture 
methylation-dependent  chromatin  silencing 
regulation  of  chromatin  assembly  or  disassembly 
ATP-dependent  chromatin  remodeling 

nucleobase,  nucleoside,  nucleotide  and  nucleic  acid  metabolic  process 

cyclic  nucleotide  biosynthetic  process 

deoxyribonucleotide  biosynthetic  process 

deoxyribonucleotide  catabolic  process 

deoxyribonucleotide  metabolic  process 

cyclic  nucleotide  catabolic  process 

negative  regulation  of  nucleotide  metabolic  process 

nucleotide  phosphorylation 
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APPENDIX  1  (CONTINUED) 


DNA  Damage,  Replication,  and  Repair  (Continued) 

nucleotide-excision  repair 

nucleotide-excision  repair,  DNA  damage  removal 
nueieotide-exeision  repair,  DNA  duplex  unwinding 
nueleotide-exeision  repair,  DNA  gap  filling 
nueieotide-exeision  repair,  DNA  ineision,  3 '-to  lesion 
nueleotide-exeision  repair,  DNA  ineision,  5 '-to  lesion 
purine  nueleotide  biosynthetie  proeess 
purine  nueleotide  metabolie  proeess 
purine  ribonueleotide  biosynthetic  process 
purine  ribonucleotide  transport 
pyridine  nueleotide  biosynthetie  proeess 
transeription-eoupled  nueleotide-exeision  repair 

Cell  Cycle 

DNA  damage  response,  signal  transduetion  by  p53  elass  mediator  resulting  in  eell  eyele 
arrest 

G1  phase  of  mitotie  eell  eyele 
G2  phase  of  mitotie  cell  eyele 
M  phase  of  mitotie  eell  eyele 
meiotie  eell  eyele 
mitotie  eell  cycle 
mitotic  cell  eyele  eheckpoint 
mitotie  eell  eyele  G2 

mitotie  eell  eyele  spindle  assembly  eheekpoint 

negative  regulation  of  cell  eyele 

negative  regulation  of  mitotie  eell  eyele 

negative  regulation  of  S  phase  of  mitotic  cell  eyele 

positive  regulation  of  eell  eyele 

re-entry  into  mitotie  eell  eyele 

regulation  of  eell  eyele 

regulation  of  mitotie  eell  eyele 

regulation  of  S  phase  of  mitotie  cell  eyele 

S  phase  of  mitotic  cell  eyele 

traversing  start  eontrol  point  of  mitotie  eell  eyele 

M  transition  of  mitotie  eell  eyele 

S  transition  of  mitotic  cell  eyele 

eell  eyele 

eell  eyele  arrest 

eell  eyele  eheekpoint 

eell  eyele  proeess 

cell  proliferation 

negative  regulation  of  eell  proliferation 
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APPENDIX  1  (CONTINUED) 


Cell  Cycle  (Continued) 

positive  regulation  of  eell  proliferation 
regulation  of  eell  proliferation 

general  transeription  from  RNA  polymerase  II  promoter 
negative  regulation  of  gene-speeifie  transeription 
negative  regulation  of  transeription 

negative  regulation  of  transcription  from  RNA  polymerase  II  promoter 
negative  regulation  of  transcription  from  RNA  polymerase  III  promoter 
negative  regulation  of  transcription,  DNA-dependent 
positive  regulation  of  gene-specific  transcription 

positive  regulation  of  specific  transcription  from  RNA  polymerase  II  promoter 
positive  regulation  of  specific  transcription  from  RNA  polymerase  II  promoter 
positive  regulation  of  transcription  from  RNA  polymerase  I  promoter 
positive  regulation  of  transcription  from  RNA  polymerase  II  promoter 
positive  regulation  of  transcription  from  RNA  polymerase  II  promoter,  global 
positive  regulation  of  transcription,  DNA-dependent 
regulation  of  transcription 

regulation  of  transcription  from  RNA  polymerase  I  promoter 

regulation  of  transcription  from  RNA  polymerase  II  promoter 

regulation  of  transcription,  DNA-dependent 

regulation  of  transcriptional  preinitiation  complex  assembly 

regulation  of  transcriptional  preinitiation  complex  assembly 

transcription 

transcription  from  RNA  polymerase  I  promoter 
transcription  from  RNA  polymerase  II  promoter 
transcription  from  RNA  polymerase  III  promoter 
transcription  initiation 

transcription  initiation  from  RNA  polymerase  I  promoter 
transcription  initiation  from  RNA  polymerase  II  promoter 
transcription  of  nuclear  rRNA  large  RNA  polymerase  I  transcript 
transcription  termination 
transcription,  DNA-dependent 
transcriptional  preinitiation  complex  assembly 

Oxidative  Stress 

heme  oxidation 

induction  of  apoptosis  by  oxidative  stress 
NADH  oxidation 

negative  regulation  of  plasma  lipoprotein  oxidation 
oxidation  reduction 
oxidative  phosphorylation 

regulation  of  transcription  from  RNA  polymerase  II  promoter  in  response  to  oxidative 
stress 
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APPENDIX  1  (CONTINUED) 


Oxidative  Stress  (Continued) 

response  to  oxidative  stress 
cell  redox  homeostasis 
redox  signal  response 
hydrogen  peroxide  catabolic  process 
removal  of  superoxide  radicals 
response  to  hydrogen  peroxide 
response  to  superoxide 
removal  of  superoxide  radicals 
response  to  oxygen  radical 
glutathione  biosynthetic  process 
glutathione  metabolic  process 
transmembrane  glutathione  transport 

Ion  Homeostasis 

response  to  zinc  ion 

sodium  ion  transport 

transition  metal  ion  transport 

zinc  ion  transport 

calcium  ion  homeostasis 

calcium  ion  transport 

cellular  calcium  ion  homeostasis 

cellular  chloride  ion  homeostasis 

cellular  iron  ion  homeostasis 

cellular  metal  ion  homeostasis 

cellular  potassium  ion  homeostasis 

cellular  response  to  iron  ion  starvation 

cellular  response  to  potassium  ion  starvation 

cellular  sodium  ion  homeostasis 

cellular  zinc  ion  homeostasis 

cobalt  ion  transport 

cytosolic  calcium  ion  homeostasis 

detection  of  calcium  ion 

elevation  of  cytosolic  calcium  ion  concentration 

endoplasmic  reticulum  calcium  ion  homeostasis 

ion  transport 

iron  ion  homeostasis 

iron  ion  transport 

magnesium  ion  transport 

manganese  ion  transport 

metal  ion  transport 

mitochondrial  calcium  ion  transport 

mitochondrial  iron  ion  transport 
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APPENDIX  1  (CONTINUED) 


Ion  Homeostasis  (Continued) 

negative  regulation  of  sodium  ion  transport 
niekel  ion  transport 

positive  regulation  of  ealeium  ion  transport 
positive  regulation  of  ealeium  ion  transport  into  eytosol 

positive  regulation  of  ealeium  ion  transport  via  store-operated  ealeium  ehannel  aetivity 

positive  regulation  of  potassium  ion  transport 

positive  regulation  of  sodium  ion  transport 

potassium  ion  import 

potassium  ion  transport 

reduetion  of  eytosolie  ealeium  ion  eoneentration 
regulation  of  sodium  ion  transport 

Apoptosis 

anti-apoptosis 

apoptosis 

apoptosis  in  response  to  endoplasmie  retieulum  stress 

DNA  damage  response,  signal  transduetion  by  p53  elass  mediator  resulting  in  induetion 
of  apoptosis 

DNA  damage  response,  signal  transduetion  resulting  in  induetion  of  apoptosis 
DNA  fragmentation  involved  in  apoptosis 
induetion  of  apoptosis 

induetion  of  apoptosis  by  extraeellular  signals 
induetion  of  apoptosis  by  intraeellular  signals 
induetion  of  apoptosis  by  oxidative  stress 
induetion  of  apoptosis  via  death  domain  reeeptors 
negative  regulation  of  anti-apoptosis 
negative  regulation  of  apoptosis 
positive  regulation  of  anti-apoptosis 
positive  regulation  of  apoptosis 

protein  insertion  into  mitoehondrial  membrane  during  induetion  of  apoptosis 

regulation  of  anti-apoptosis 

regulation  of  apoptosis 

transformed  eell  apoptosis 

aetivation  of  pro-apoptotie  gene  produets 

eell  death 

death 

positive  regulation  of  programmed  eell  death 
programmed  eell  death 

Energy  regulation 

eellular  glueose  homeostasis 
glueose  1 -phosphate  metabolie  proeess 
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Energy  Regulation  (Continued) 

glucose  6-phosphate  metabolie  proeess 
glueose  eatabolie  proeess 
glueose  homeostasis 
glueose  metabolie  proeess 
glueose  transport 
gluoose-6-phosphate  transport 
negative  regulation  of  glueose  import 
positive  regulation  of  glueose  import 
glyeolysis 

regulation  of  glyeolysis 

Hypoxic  response 

response  to  hypoxia 
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APPENDIX  2:  PROBE  SETS  DIFFERENTIALLY  EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1367576_at 

Gpxl 

glutathione  peroxidase  1 

Blue 

1367586_at 

Ldha 

laetate  dehydrogenase  A 

Blue 

1367603_at 

Tpil 

triosephosphate  isomerase  1 

Blue 

1367705_at 

Glrxl 

glutaredoxin  1  (thioltransferase) 

Yellow 

1367707_at 

Fasn 

fatty  aeid  synthase 

Yellow 

1367708_a_at 

Fasn 

fatty  aeid  synthase 

Yellow 

1367722_at 

Dpp7 

dipeptidylpeptidase  7 

Yellow 

1367754_s_at 

Asl 

argininosuecinate  lyase 

Yellow 

1367758_at 

Afp 

alpha- fetoprotein 

Yellow 

1367760_at 

Map2kl 

mitogen  activated  protein  kinase  kinase  1 

Blue 

1367775_at 

Amacr 

alpha-methylacyl-CoA  racemase 

Yellow 

1367786_at 

Psmb8 

proteasome  (prosome,  macropain)  subunit, 
beta  type  8 

Yellow 

1367866_at 

Fbln5 

fibulin  5 

Yellow 

1367898_at 

Bnip31 

BCL2/adenovirus  ElB  interacting  protein 

3 -like 

Blue 

1 36790  l_at 

Gusb 

glucuronidase,  beta 

Yellow 

1368000_at 

C3 

complement  component  3 

Blue 

1368025_at 

Ddit4 

DNA-damage-inducible  transcript  4 

Blue 

1368079_at 

Pdkl 

pyruvate  dehydrogenase  kinase,  isozyme  1 

Blue 

1368082_at 

Slc4a2 

solute  carrier  family  4  (anion  exchanger), 
member  2 

Yellow 

1368130_at 

Aldh3al 

aldehyde  dehydrogenase  3  family,  member 
A1 

Blue 

1368174_at 

Egln3 

EGE  nine  homolog  3  (C.  elegans) 

Blue 

1 36824  l_a_at 

Flotl 

flotillin  1 

Yellow 

1368258_at 

Apln 

apelin 

Blue 

1368280_at 

Ctsc 

cathepsin  C 

Yellow 

1368308_at 

Myc 

myelocytomatosis  oncogene 

Blue 

1368428_at 

Xpnpep2 

X-prolyl  aminopeptidase  (aminopeptidase 

P)  2,  membrane-bound 

Yellow 

1368488_at 

Nfd3 

nuclear  factor,  interleukin  3  regulated 

Blue 

136851  l_at 

Bhlhe41 

basic  helix-loop-helix  family,  member  e41 

Blue 

13685 19_at 

Serpinel 

serine  (or  cysteine)  peptidase  inhibitor, 
clade  E,  member  1 

Blue 

1368618_at 

Grbl4 

growth  factor  receptor  bound  protein  14 

Yellow 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1368674_at 

Pygl 

pbospborylase,  glycogen,  liver 

Blue 

1368835_at 

Stall 

signal  transducer  and  activator  of 
transcription  1 

Yellow 

1368863_at 

Nme3 

non-metastatic  cells  3,  protein  expressed  in 

Yellow 

1368916_at 

Asl 

argininosuccinate  lyase 

Yellow 

1369224_at 

Cdhl7 

cadberin  17 

Yellow 

1369473_at 

Pgml 

pbospboglucomutase  1 

Blue 

1369788_s_at 

Jun 

Jun  oncogene 

Blue 

1369958_at 

Rhob 

ras  bomolog  gene  family,  member  B 

Blue 

1370030_at 

Gclm 

glutamate  cysteine  ligase,  modifier  subunit 

Blue 

1370062_at 

Higdla 

HlGl  domain  family,  member  lA 

Blue 

1370120_at 

Fstl3 

follistatin-like  3  (secreted  glycoprotein) 

Blue 

1370138_at 

Lefl 

lympboid  enhancer  binding  factor  1 

Yellow 

1370235_at 

Dbi 

diazepam  binding  inhibitor  (GABA 
receptor  modulator,  acyl-Coenzyme  A 
binding  pro 

Yellow 

1370237_at 

Hadb 

hydroxyacyl-Coenzyme  A  dehydrogenase 

Yellow 

1370688_at 

Gclc 

glutamate-cysteine  ligase,  catalytic  subunit 

Blue 

1 37082  l_at 

Tpmt 

thiopurine  S-methyltransferase 

Yellow 

1370848_at 

Slc2al 

solute  earrier  family  2  (facilitated  glucose 
transporter),  member  1 

Blue 

1370954_at 

P4bal 

procollagen-proline,  2-oxoglutarate  4- 
dioxygenase  (proline  4-hydroxylase), 
alpha 

Blue 

1370959_at 

Col3al 

collagen,  type  111,  alpha  1 

Yellow 

1370975_at 

Jmjdla 

jumonji  domain  containing  lA 

Blue 

1371005_at 

Abed 

ATP -binding  cassette,  sub-family  C 
(CFTR/MRP),  member  1 

Blue 

1371467_at 

LOC293103 

similar  to  RIKEN  cDNA  0610007P06 

Blue 

1371824_at 

Ak311 

adenylate  kinase  3 -like  1 

Blue 

1372012_at 

Dbcr24 

24-dehydrocholesterol  reductase 

Yellow 

1372178_at 

— 

— 

Green 

1 37226  l_at 

— 

— 

Blue 

1372277_at 

LOC641316 

similar  to  aldehyde  dehydrogenase  4 
family,  member  A1 

Yellow 

1372318  at 

— 

— 

Yellow 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1372602_at 

Stbdl 

starch  binding  domain  1 

Blue 

1372610_at 

P4ha2 

procollagen-proline,  2-oxoglutarate  4- 
dioxygenase  (proline  4-hydroxylase) 

Blue 

1372653_at 

Fkbpll 

FK506  binding  protein  1 1 

Blue 

1372706_at 

Hexb 

hexosaminidase  B 

Blue 

1372760_at 

— 

— 

Blue 

1372830_at 

Arlbipl 

ADP-ribosylation  faetor-like  6  interaeting 
protein  1 

Blue 

1372863_at 

Mycbp2 

MYC  binding  protein  2 

Blue 

1372993_at 

— 

— 

Blue 

1373026_at 

Spc24 

SPC24,  NDC80  kinetoehore  eomplex 
eomponent,  homolog  (S.  eerevisiae) 

White 

1373036_at 

RGD1561455 

similar  to  Ras  GTPase-aetivating-like 
protein  IQGAP2 

Yellow 

1373037_at 

Ube216 

ubiquitin-conjugating  enzyme  E2L  6 

Yellow 

1373093_at 

Errfil 

ERBB  reeeptor  feedbaek  inhibitor  1 

Blue 

1373108_at 

Ppplr3c 

protein  phosphatase  1 ,  regulatory 
(inhibitor)  subunit  3C 

Blue 

1373433_at 

Nsbpl 

nueleosomal  binding  protein  1 

Yellow 

1373438_at 

Ube2o 

ubiquitin-eonjugating  enzyme  E20 

Blue 

1373475_at 

Ccdc58 

eoiled-coil  domain  eontaining  58 

Blue 

1373524_at 

... 

— 

Blue 

1373544_at 

Cxcl9 

ehemokine  (C-X-C  motif)  ligand  9 

Yellow 

1373685_at 

Ankrd37 

ankyrin  repeat  domain  37 

Blue 

1373736_at 

— 

— 

Blue 

1373742_at 

Spsb2 

splA/ryanodine  reeeptor  domain  and  SOCS 
box  eontaining  2 

Yellow 

1373807_at 

Vegfa 

vaseular  endothelial  growth  factor  A 

Blue 

1373860_at 

Sox4 

SRY  (sex  determining  region  Y)-box  4 

Blue 

1373942_at 

... 

— 

White 

1374105_at 

Higdla 

HIGl  domain  family,  member  lA 

Blue 

1374244_at 

— 

— 

Yellow 

1374433_at 

... 

— 

Blue 

1374447_at 

Usp9x 

ubiquitin  speeific  peptidase  9,  X-linked 

Blue 

1374524_at 

Scly 

selenocysteine  lyase 

Yellow 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1374575_at 

Creb311 

cAMP  responsive  element  binding  protein 

3 -like  1 

Yellow 

1374632_at 

Jmjd6 

jumonji  domain  eontaining  6 

Blue 

1374663_at 

— 

— 

Blue 

1374778_at 

Ctsc 

cathepsin  C 

Yellow 

1374828_at 

Pdia5 

protein  disulfide  isomerase  family  A, 
member  5 

Blue 

1374883_at 

Mtmr7 

myotubularin  related  protein  7 

Yellow 

1374906_at 

Rnfl  1 3a1 

ring  finger  protein  1 13A1 

Blue 

1374932_at 

— 

— 

Blue 

1374948_at 

Tmeml06a 

transmembrane  protein  106A 

Blue 

1375123_at 

Sox4 

SRY  (sex  determining  region  Y)-box  4 

Blue 

1375230_at 

— 

— 

Blue 

1375247_at 

Mgll 

Monoglyceride  lipase 

Yellow 

1375362_at 

Sppl2a 

signal  peptide  peptidase-like  2A 

Yellow 

1375374_at 

— 

— 

Green 

1375388_at 

— 

— 

Blue 

1375532_at 

— 

— 

Yellow 

1375793_at 

Mill 

myeloid/lymphoid  or  mixed-lineage 
leukemia  1 

Blue 

1375796_at 

Igtp 

interferon  gamma  induced  GTPase 

Yellow 

137615 l_a_at 

... 

— 

Yellow 

1376163_at 

Clrn3 

clarin  3 

Yellow 

1376677_a_at 

— 

— 

Blue 

1376724_at 

... 

— 

Blue 

1376917_at 

Znf292 

zinc  finger  protein  292 

Blue 

1377042_at 

Pcgf5 

polycomb  group  ring  finger  5 

Blue 

1377043_at 

— 

— 

Blue 

1377234_at 

— 

— 

Blue 

1377277_at 

... 

— 

Yellow 

137755  l_at 

... 

— 

Blue 

1377705_at 

... 

— 

Blue 

1377706_x_at 

— 

— 

Blue 

1377818  at 

— 

— 

Blue 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1377928_at 

Tm4sf20 

transmembrane  4  L  six  family  member  20 

Yellow 

1 37808  l_at 

... 

— 

Blue 

1378144_at 

Kank4 

KN  motif  and  ankyrin  repeat  domains  4 

Yellow 

1378210_at 

Mill 

myeloid/lymphoid  or  mixed-lineage 
leukemia  1 

Blue 

1 37826  l_at 

— 

— 

Blue 

13783 14_at 

— 

— 

Blue 

1378361_at 

Chd7 

chromodomain  helicase  DNA  binding 
protein  7 

Blue 

1378413_at 

— 

— 

Yellow 

1378479_at 

Cln3 

ceroid-lipofuscinosis,  neuronal  3 

Yellow 

1378808_at 

— 

— 

Blue 

1378833_at 

— 

— 

Blue 

1379356_at 

... 

— 

Yellow 

1379358_at 

Samd4a 

Sterile  alpha  motif  domain  containing  4A 

Blue 

1379365_at 

Cxclll 

chemokine  (C-X-C  motif)  ligand  1 1 

Yellow 

1379415_at 

— 

— 

Blue 

1379440_at 

Fstl3 

follistatin-like  3  (secreted  glycoprotein) 

Blue 

1379452_at 

Gas2 

growth  arrest-specific  2 

Yellow 

1379498_at 

... 

— 

Blue 

1379580_at 

Clid6 

chromodomain  helicase  DNA  binding 
protein  6 

Blue 

1379636_at 

Fam82a 

family  with  sequence  similarity  82, 
member  A 

Yellow 

1379739_at 

— 

— 

Blue 

1379900_at 

... 

— 

Yellow 

1379932_at 

Clcn4-2 

chloride  channel  4-2 

Yellow 

1379990_at 

... 

— 

Blue 

1380474_at 

Loxl2 

Lysyl  oxidase-like  2 

Blue 

1380794_at 

— 

— 

Blue 

1380997_at 

Ccnj 

cyclin  J 

Blue 

1381229_at 

— 

— 

Blue 

1381684_at 

... 

— 

Blue 

1381823_at 

— 

— 

Blue 

1381852  at 

— 

— 

Yellow 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1382001_at 

LOC502710 

similar  to  Myeloid/lymphoid  or  mixed- 
lineage  leukemia  protein  3  homolog 

Blue 

1382059_at 

Fbxo30 

F-box  protein  30 

Blue 

1382072_at 

01fml2a 

olfactomedin-like  2A 

Yellow 

1382255_at 

... 

— 

Blue 

1382266_at 

Gprl46 

G  protein-coupled  receptor  146 

Blue 

1382325_at 

Gcat 

glycine  C-acetyltransferase  (2-ammo-3- 
ketobutyrate-coenzyme  A  ligase) 

Yellow 

1382348_at 

Dhodh 

dihydroorotate  dehydrogenase 

Blue 

1382400_at 

Rif 

rearranged  L-myc  fusion 

Blue 

1382658_at 

LOC502710 

similar  to  Myeloid/lymphoid  or  mixed- 
lineage  leukemia  protein  3  homolog 
(Histone 

Blue 

1382776_at 

Mboatl 

membrane  bound  0-acyltransferase 
domain  containing  1 

Yellow 

1382797_at 

RGD 1560433 

similar  to  1 5000 19C06Rik  protein 

Blue 

1382924_at 

Pankl 

pantothenate  kinase  1 

Yellow 

1382932_at 

Lrrc50 

leucine  rich  repeat  containing  50 

Blue 

1382942_at 

... 

— 

Blue 

1382982_at 

— 

— 

Blue 

1383013_at 

Klfl3 

Kruppel-like  factor  13 

Blue 

1383058_at 

— 

— 

Green 

1383401_at 

Tes 

testis  derived  transcript 

Blue 

1383409_at 

... 

— 

Blue 

1383432_at 

Dapll 

death  associated  protein-like  1 

Yellow 

1383500_at 

Rrbpl 

ribosome  binding  protein  1 

Blue 

1383785_at 

Lefl 

Lymphoid  enhancer  binding  factor  1 

Yellow 

1383826_at 

Rab40b 

Rab40b,  member  RAS  oncogene  family 

Blue 

1383939_at 

... 

— 

Blue 

1384336_at 

Cln3 

ceroid-lipofuscinosis,  neuronal  3 

Yellow 

1384406_at 

LOC685243 

similar  to  C-C  chemokine  receptor  type  1 1 
(C-C  CKR-1 1)  (CC-CKR-1 1)  (CCR-1 1) 

Yellow 

1384432_at 

— 

— 

Blue 

1384475_at 

... 

— 

Yellow 

1384971_at 

Depdc7 

DEP  domain  containing  7 

Yellow 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1385050_at 

Znf292 

zinc  finger  protein  292 

Blue 

1385166_at 

... 

— 

Yellow 

1385171_s_at 

Smtnl2 

smoothelin-like  2 

Blue 

1385334_at 

... 

— 

Yellow 

1385706_at 

— 

— 

Blue 

1385827_at 

Clcfl 

eardiotrophin-like  eytokine  faetor  1 

Green 

138623  l_at 

RGD 1306474 

similar  to  RIKEN  eDNA  9530003J23 

Yellow 

1386352_at 

... 

— 

Blue 

1386695_at 

... 

— 

Yellow 

1386716_at 

— 

— 

Yellow 

1386773_at 

Bid 

biotinidase 

Yellow 

1386893_at 

Grn 

granulin 

Yellow 

1386908_at 

Glrxl 

glutaredoxin  1  (thioltransferase) 

Yellow 

1386946_at 

Cptla 

earnitine  palmitoyltransferase  la,  liver 

Yellow 

1386958_at 

Txnrdl 

thioredoxin  reductase  1 

Blue 

1386978_at 

Bnip31 

BCL2/adenovirus  ElB  interaeting  protein 

3 -like 

Blue 

1387028_a_at 

Idl 

inhibitor  of  DNA  binding  1 

Blue 

1387058_at 

Pctp 

phosphatidyleholine  transfer  protein 

Yellow 

1387060_at 

Klf6 

Kruppel-like  faetor  6 

Green 

1387121_a_at 

Ndrg2 

NDRG  family  member  2 

Blue 

1387130_at 

Slc40al 

solute  earrier  family  39  (iron-regulated 
transporter),  member  1 

Blue 

1387145_at 

Gjbl 

gap  junetion  protein,  beta  1 

Yellow 

1387178_a_at 

Cbs 

eystathionine  beta  synthase 

Yellow 

1387219_at 

Adm 

adrenomedullin 

Blue 

1387271_at 

Phyh 

phytanoyl-CoA  2-hydroxylase 

Yellow 

13873 14_at 

Sultlbl 

sulfotransferase  family,  eytosolie,  IB, 
member  1 

Yellow 

1387338_s_at 

Bcl2111 

BCE2-like  1 1  (apoptosis  faeilitator) 

Blue 

1387354_at 

Stall 

signal  transdueer  and  aetivator  of 
transeription  1 

Yellow 

1387361_s_at 

Pgkl 

phosphoglyeerate  kinase  1 

Blue 

1387528_at 

Mbl2 

mannose-binding  lectin  (protein  C)  2 

Yellow 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1387805_at 

Bnip3 

BCE2/adenovirus  ElB  19  kDa- interacting 
protein  3 

Blue 

1387893_at 

Cls 

eomplement  eomponent  1 ,  s  subeomponent 

Y  ellow 

1387914_at 

Cyp27al 

cytoehrome  P450,  family  27,  subfamily  a, 
polypeptide  1 

Yellow 

1387964_a_at 

Eroll 

EROl-like  (S.  eerevisiae) 

Blue 

1387995_a_at 

Ilitm3 

interferon  indueed  transmembrane  protein 

3 

Yellow 

1388102_at 

Ptgrl 

prostaglandin  reduetase  1 

Blue 

1388108_at 

Elovl6 

EEOVE  family  member  6,  elongation  of 
long  ehain  fatty  aeids  (yeast) 

Yellow 

1388318_at 

Pgkl 

phosphoglyeerate  kinase  1 

Blue 

1388587_at 

Ier3 

immediate  early  response  3 

Blue 

1388600_at 

— 

— 

Yellow 

1388634_at 

Pgml 

phosphoglucomutase  1 

Blue 

1388666_at 

Enel 

eetodermal-neural  eortex  1 

Blue 

1388742_at 

... 

— 

Blue 

1388872_at 

... 

— 

Yellow 

1388948_at 

StardlO 

StAR-related  lipid  transfer  (START) 
domain  containing  10 

Yellow 

1388985_at 

— 

— 

Blue 

1388986_at 

... 

— 

Green 

1389014_at 

Nampt 

nicotinamide  phosphoribosyltransferase 

Blue 

1389075_at 

— 

— 

Blue 

1389132_at 

Hipl 

huntingtin  interaeting  protein  1 

Yellow 

1389199_at 

RGD 1309079 

similar  to  Ab2-095 

Blue 

1389207_at 

Eglnl 

EGL  nine  homolog  1  (C.  elegans) 

Blue 

1389230_at 

— 

— 

Blue 

1389297_at 

Eroll 

EROl-like  (S.  eerevisiae) 

Blue 

13893 16_at 

Usp9x 

ubiquitin  speeifie  peptidase  9,  X-linked 

Blue 

1389409_at 

Tes 

testis  derived  transeript 

Blue 

138953  l_at 

Zfp330 

zinc  finger  protein  330 

Green 

1389555_at 

Tcfl9 

transeription  factor  19 

White 

1389618  at 

— 

— 

Blue 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

Gene 

Gene  Title 

Cluster 

ID 

Symbol 

1389797_at 

— 

— 

Blue 

1389858_at 

Tkl 

thymidine  kinase  1 ,  soluble 

White 

1390020_at 

Ogdh 

oxoglutarate  (alpha-ketoglutarate) 

Yellow 

dehydrogenase  (lipoamide) 

1390082_at 

Hipl 

huntingtin  interaeting  protein  1 

Yellow 

1390109_at 

— 

— 

Blue 

1390203_at 

— 

— 

Blue 

1390394_at 

LOC500034 

similar  to  CG3570-PA 

Blue 

1390412_at 

Slc40al 

solute  earrier  family  39  (iron-regulated 

Blue 

transporter),  member  1 

1390423_at 

Mycbp2 

MYC  binding  protein  2 

Blue 

1390530_at 

— 

— 

Yellow 

1390532_at 

— 

— 

Yellow 

1 39059  l_at 

Slcl7a3 

solute  earrier  family  17  (sodium 

Blue 

phosphate),  member  3 

1390709_at 

Trio 

triple  funetional  domain  (PTPRF 

Blue 

interaeting) 

1390723_at 

— 

— 

Blue 

1390862_at 

Znf629 

zine  finger  protein  629 

Blue 

1390932_at 

— 

— 

Blue 

1390945_at 

Znf292 

zine  finger  protein  292 

Blue 

1390977_at 

St71 

suppression  of  tumorigenieity  7-like 

Blue 

1390993_at 

Pbld 

phenazine  biosynthesis-like  protein  domain 

Yellow 

eontaining 

1391212_at 

Tceall 

transeription  elongation  faetor  A  (Sll)-like 

1 

Yellow 

1391215_at 

— 

— 

Blue 

1391323_at 

Srprb  -  Tf 

signal  reeognition  partiele  reeeptor,  B 

Blue 

subunit  -  transferrin 

1391458_at 

Ndrgl 

N-mye  downstream  regulated  gene  1 

Blue 

1391675_at 

RGD 1309708 

similar  to  RIKEN  oDNA  4930455F23 

Blue 

1391683_at 

— 

— 

Blue 

1391880_at 

— 

— 

Blue 

1392038_at 

Arlbipl 

ADP-ribosylation  faetor-like  6  interaeting 

Blue 

protein  1 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1392603_at 

— 

— 

Blue 

1392633_at 

... 

— 

Blue 

1392736_at 

... 

— 

Blue 

1392747_at 

... 

— 

Blue 

1 39294  l_at 

Hipl 

huntingtin  interacting  protein  1 

Yellow 

1392946_at 

— 

— 

Yellow 

1392952_at 

Acsf2 

acyl-CoA  synthetase  family  member  2 

Yellow 

1393047_at 

Rnf217 

ring  linger  protein  217 

Blue 

1393236_at 

Riok3 

RIO  kinase  3  (yeast) 

Blue 

1393245_at 

Phyh 

Phytanoyl-CoA  2-hydroxylase 

Yellow 

1393377_at 

— 

— 

Yellow 

1393492_at 

... 

— 

Blue 

13935 16_at 

... 

— 

Yellow 

1393612_a_at 

Depdc7 

DEP  domain  containing  7 

Yellow 

1393627_at 

Osap 

ovary-specific  acidic  protein 

Blue 

1 39440  l_at 

Elovl6 

EEOVE  family  member  6,  elongation  of 
long  chain  fatty  acids  (yeast) 

Yellow 

1394459_at 

FAM120C 

family  with  sequence  similarity  120C 

Yellow 

1394844_s_at 

Cyp4a2  - 
Cyp4a3 

cytochrome  P450,  family  4,  subfamily  a, 
polypeptide  2  -  cytochrome  P450,  family 

Yellow 

1394948_at 

— 

— 

Blue 

1 39504  l_at 

Mttp 

microsomal  triglyceride  transfer  protein 

Yellow 

1395236_at 

Ppplr3c 

protein  phosphatase  1 ,  regulatory 
(inhibitor)  subunit  3C 

Blue 

1395427_at 

— 

— 

White 

13955 19_at 

MGC95152 

similar  to  B230212E03Rik  protein 

Blue 

1395744_at 

... 

— 

Yellow 

1395896_at 

— 

— 

Blue 

1396262_at 

Nampt 

nicotinamide  phosphoribosyltransferase 

Blue 

1396539_at 

— 

— 

Blue 

1397304_at 

Igtp 

interferon  gamma  induced  GTPase 

Yellow 

13973  lOat 

— 

— 

Blue 

1397443_at 

— 

— 

Blue 

1397633_at 

— 

— 

Green 
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APPENDIX  2  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL 


Affymetrix  Gene  Gene  Title  Cluster 


ID _ Symbol 


1398380_at 

Vwal 

von  Willebrand  factor  A  domain 
containing  1 

Green 

139841  l_at 

— 

— 

Blue 

1398727_at 

— 

— 

Yellow 

1 39879  l_at 

Txnrdl 

thioredoxin  reductase  1 

Blue 

1398923_at 

Rnasek 

ribonuelease,  RNase  K 

Yellow 

13991 19_at 

Znf292 

zine  finger  protein  292 

Blue 

1397838_at 

RGD 1564792 

RGD  1564792 

Blue 

1 37454  l_at 

Gysl 

glyeogen  synthase  1,  muscle 

Blue 

1387671_at 

Sctr 

secretin  reeeptor 

Yellow 

138965  l_at 

Apln 

apelin 

Blue 

1391958_at 

— 

— 

Blue 

1372542_at 

... 

— 

Blue 

1372523_at 

Gclc 

glutamate-eysteine  ligase,  catalytie  subunit 

Blue 

1374537_at 

Chsyl 

chondroitin  sulfate  synthase  1 

Blue 

1377902_a_at 

Rad52 

RAD52  homolog  (S.  cerevisiae) 

Blue 

1381902_at 

Znf292 

zinc  finger  protein  292 

Blue 

1385017_at 

Aqr  -  Znf770 

aquarius  homolog  (mouse)  -  zinc  finger 

Blue 

1387622_at 

Rfng 

RFNG  0-fueosylpeptide  3-beta-N- 
acetylglucosaminyltransferase 

Yellow 

1389641_at 

Snapc4 

small  nuclear  RNA  aetivating  eomplex, 
polypeptide  4 

Blue 

1392264_s_at 

Serpinel 

serine  (or  cysteine)  peptidase  inhibitor, 
clade  E,  member  1 

Blue 

1397709_at 

Arl6ipl 

ADP-ribosylation  faetor-like  6  interaeting 
protein  1 

Blue 

74 


APPENDIX  3:  PROBE  SETS  DIFFERENTIALLY  EXPRESSED  BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1 36767  l_at 

Pena 

proliferating  cell  nuclear  antigen 

White 

1367790_at 

Sndl 

staphylocoecal  nuclease  and  tudor  domain 
containing  1 

Yellow 

1367794_at 

A2m 

alpha-2-macroglobulin 

Yellow 

1367847_at 

Nuprl 

nuclear  protein  1 

Yellow 

1367917_at 

Cyp2d2 

eytochrome  P450,  family  2,  subfamily  d, 
polypeptide  2 

Yellow 

1367983_at 

Fenl 

flap  strueture-speeific  endonuclease  1 

White 

1367994_at 

Dpyd 

dihydropyrimidine  dehydrogenase 

Yellow 

1368060_at 

Hrspl2 

heat-responsive  protein  12 

Yellow 

1368092_at 

Fah 

fumarylaeetoacetate  hydrolase 

Yellow 

1368204_at 

Ligl 

ligase  I,  DNA,  ATP-dependent 

White 

1368224_at 

Serpina3n 

serine  (or  cysteine)  peptidase  inhibitor,  elade 
A,  member  3N 

Yellow 

1368244_at 

As3mt 

arsenic  (+3  oxidation  state) 
methyltransferase 

Yellow 

1368253_at 

Gamt 

guanidinoaeetate  N-methyltransferase 

Yellow 

1368387_at 

Bdhl 

3-hydroxybutyrate  dehydrogenase,  type  1 

Yellow 

1368399_a_at 

Pgcp 

plasma  glutamate  earboxypeptidase 

Yellow 

1368452_at 

Abcc6 

ATP -binding  eassette,  sub-family  C 
(CFTR/MRP),  member  6 

Yellow 

1368497_at 

Abcc2 

ATP -binding  eassette,  sub-family  C 
(CFTR/MRP),  member  2 

Green 

1368569_at 

Akrlb7 

aldo-keto  reductase  family  1 ,  member  B7 

Yellow 

1368642_at 

Cdh2 

cadherin  2 

Yellow 

1368658_at 

Cntf 

eiliary  neurotrophic  factor 

White 

1368707_at 

Itih4 

inter  alpha-trypsin  inhibitor,  heavy  chain  4 

Yellow 

1368860_at 

Phldal 

pleckstrin  homology-like  domain,  family  A, 
member  1 

Yellow 

1368924_at 

Ghr 

growth  hormone  receptor 

Yellow 

1369086  a  at 

Cacnald 

calcium  channel,  voltage-dependent,  L  type, 
alpha  ID  subunit 

Yellow 

13691 7  l_at 

Mstl 

Maerophage  stimulating  1  (hepatoeyte 
growth  faetor-like) 

Yellow 

1369225_at 

Kng2 

kininogen  2 

Yellow 

13693 18_at 

Fhit 

fragile  histidine  triad  gene 

Yellow 

1369610_at 

Lin7c 

lin-7  homolog  C  (C.  elegans) 

White 

1369629  at 

Adk 

adenosine  kinase 

Yellow 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1369777_a_at 

Shank2 

SH3  and  multiple  ankyrin  repeat  domains  2 

Yellow 

1369837_at 

Gulo 

gulonolaetone  (L-)  oxidase 

Yellow 

1369852_at 

FIO 

eoagulation  faetor  X 

Yellow 

1369909_s_at 

Tmeml50 

transmembrane  protein  150 

Yellow 

1369940_at 

Taldol 

transaldolase  1 

White 

1370127_at 

Poldl 

polymerase  (DNA  directed),  delta  1, 
catalytic  subunit 

White 

1370336_at 

Osginl 

oxidative  stress  induced  growth  inhibitor  1 

Blue 

1370352_at 

Cesll  ///  Es22 
///rCG_44273 

carboxylesterase-like  1  III  esterase  22 
///carboxylesterase  ES-4 

White 

1370436_at 

Acsm2 

acyl-CoA  synthetase  medium-chain  family 
member  2 

Yellow 

1370475_at 

Cyp2b3 

cytochrome  P450IIB3 

Yellow 

1370599_a_at 

Ptprd 

protein  tyrosine  phosphatase,  receptor  type, 

D 

UDP-glucuronosyltransferase 

Yellow 

1370615_at 

LOC286989 

White 

1370693_a_at 

Cnp 

2',3'-cyclic  nucleotide  3'  phosphodiesterase 

Yellow 

1370979_at 

Ddx20 

DEAD  (Asp-Glu- Ala-Asp)  box  polypeptide 
20 

fibrinogen  alpha  chain 

White 

1370992_a_at 

Fga 

Yellow 

1371014_at 

Plcbl 

phospholipase  C,  beta  1  (phosphoinositide- 
specific) 

Yellow 

1371069_at 

Slc22a23 

solute  carrier  family  22,  member  23 

Yellow 

1371074_a_at 

Mcm6 

minichromosome  maintenance  complex 
component  6 

White 

1371081_at 

Rapgef4 

Rap  guanine  nucleotide  exchange  factor 
(GEE)  4 

Yellow 

1371094_at 

Lhx2 

EIM  homeobox  2 

White 

1371099_at 

— 

— 

Yellow 

1371104_at 

Srebfl 

sterol  regulatory  element  binding 
transcription  factor  1 

Yellow 

1371124_a_at 

Kng2 

kininogen  2 

Yellow 

1371332_at 

LOC684681 

similar  to  Histone  HI. 2  (HI  VAR.l)  (Hie) 

Yellow 

1371541_at 

Mylk 

myosin  light  chain  kinase 

Yellow 

1371571_at 

App 

amyloid  beta  (A4)  precursor  protein 

Yellow 

1371572_at 

App 

amyloid  beta  (A4)  precursor  protein 

Yellow 

1371813_at 

Hirip3 

HIRA  interacting  protein  3 

White 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1371942_at 

Gstt3 

glutathione  S-transferase,  theta  3 

Yellow 

1371963_at 

Pcca 

propionyl-coenzyme  A  carboxylase,  alpha 
polypeptide 

Yellow 

1371970_at 

Faml  11a 

family  with  sequence  similarity  111, 
member  A 

White 

1371988_at 

Manlal 

mannosidase,  alpha,  class  lA,  member  1 

Yellow 

1371990_at 

Mpnd 

MPN  domain  containing 

Y  ellow 

1371995_at 

Klhl21 

kelch-like21  (Drosophila) 

White 

1372378_at 

Dis312 

D1S3  mitotic  control  homolog  (S. 
cerevisiae)-like  2 

Yellow 

1372406_at 

Mcm3 

minichromosome  maintenance  complex 
component  3 

White 

1372438_at 

Nit2 

nitrilase  family,  member  2 

Yellow 

13725 16_at 

Kif22 

kinesin  family  member  22 

White 

1372903_at 

Kifl8b 

kinesin  family  member  1 8B 

White 

1373079_at 

— 

— 

Yellow 

1373178_at 

— 

— 

Yellow 

1373283_at 

— 

— 

White 

1373290_at 

Ezh2 

enhancer  of  zeste  homolog  2  (Drosophila) 

White 

1373419_at 

Ptprg 

Protein  tyrosine  phosphatase,  receptor  type. 

Yellow 

1373489_at 

Chafla 

chromatin  assembly  factor  1 ,  subunit  A 
(pl50) 

White 

1373530_at 

Ccnel 

cyclin  El 

White 

1373557_at 

Mcm4 

minichromosome  maintenance  complex 
component  4 

White 

1373628_at 

— 

— 

Yellow 

1373654_at 

Anxa8 

annexin  A8 

Yellow 

1373773_at 

Gpm6a 

glycoprotein  m6a 

Yellow 

1373892_at 

— 

— 

Yellow 

1373908_at 

— 

— 

Yellow 

1373935_at 

Pold2 

polymerase  (DNA  directed),  delta  2, 
regulatory  subunit 

White 

1374036_at 

Mcm2 

minichromosome  maintenance  complex 
component  2 

White 

1374076_at 

LOC690743 
///  Rfwd3 

similar  to  mixed  lineage  kinase  domain-like 
///  ring  finger  and  WD  repeat  domain 

White 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1374222_at 

Slc22al8 

solute  carrier  family  22,  member  18 

Yellow 

1374320_at 

F5 

coagulation  factor  V  (proaccelerin,  labile 
factor) 

Yellow 

1374359_at 

Ccne2 

cyclin  E2 

White 

1374417_at 

Nrbp2 

nuclear  receptor  binding  protein  2 

Yellow 

1374432_at 

— 

— 

Yellow 

1374495_at 

Lrba 

EPS-responsive  vesicle  trafficking,  beach 
and  anchor  containing 

Yellow 

13745 12_at 

LOC684314 

similar  to  CG9886-like 

Yellow 

1374532_at 

Ptges2 

prostaglandin  E  synthase  2 

White 

1374540_at 

Cdca7 

cell  division  cycle  associated  7 

White 

1374747_at 

Pftkl 

PET  AIRE  protein  kinase  1 

Yellow 

1 37486  l_at 

Tle2 

transducin-like  enhancer  of  split  2  (E(spl) 
homolog.  Drosophila) 

Yellow 

137491 l_at 

RGD1303142 

oxidative  stress  responsive  gene 

Green 

1374915_at 

— 

— 

Yellow 

1374995_at 

Elmo3 

engulfment  and  cell  motility  3 

Yellow 

1375026_at 

Calml4 

calmodulin-like  4 

Yellow 

1375224_at 

Phlda3 

pleckstrin  homology-like  domain,  family  A, 
member  3 

White 

13753 14_at 

Bcas3 

breast  carcinoma  amplified  sequence  3 

Yellow 

1375343_at 

— 

— 

Yellow 

13755 17_at 

Trp53inp2 

tumor  protein  p53  inducible  nuclear  protein 

9 

Yellow 

1375562_at 

— 

L 

Yellow 

1375612_at 

— 

— 

Green 

1375856_at 

— 

— 

Yellow 

1375857_at 

Myof 

myoferlin 

Yellow 

1375943_at 

CpllO 

CPI  10  protein 

White 

1376018_at 

Lmnb2 

Eamin  B2 

Green 

1376026_at 

Donson 

downstream  neighbor  of  SON 

White 

1376055_at 

Mcm5 

minichromosome  maintenance  complex 
component  5 

White 

1376071_at 

— 

— 

Yellow 

1376124_at 

Iqub 

IQ  motif  and  ubiquitin  domain  containing 

White 

1376132_at 

Gtdcl 

glycosyltransferase-like  domain  containing  1 

Y  ellow 

13763 17_at 

Orc61 

origin  recognition  complex,  subunit  6  like 

White 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1376659_at 

Nubpl 

nucleotide  binding  protein-like 

Yellow 

1376687_at 

Uspl 

ubiquitin  speeific  peptidase  1 

White 

137673 l_at 

— 

— 

Yellow 

1376733_at 

Igsfll 

immunoglobulin  superfamily,  member  1 1 

Yellow 

1376747_at 

— 

— 

Yellow 

1376774_at 

Exoc4 

Exoeyst  eomplex  eomponent  4 

Yellow 

1376810_at 

Mphosph8  /// 
Parp4 

M-phase  phosphoprotein  8  III  poly  (ADP- 
ribose)  polymerase  family,  member  4 

Yellow 

1376840_at 

— 

— 

White 

1376845_at 

isg  12(b) 

putative  ISG  12(b)  protein 

Yellow 

1376860_at 

Slc22a23 

solute  earner  family  22,  member  23 

Yellow 

13771 14_at 

— 

— 

Green 

1377183_at 

— 

— 

Yellow 

1377254_a_at 

Cohhl 

Cohen  syndrome  homolog  1 

Yellow 

1377410_at 

E2f8 

E2P  transeription  faetor  8 

White 

1377496_at 

RGD 1305441 

similar  to  RIKEN  oDNA  5730453116 

Green 

1377718_at 

Eancb 

Eaneoni  anemia,  eomplementation  group  B 

White 

1377832_at 

Plk4 

polo-like  kinase  4  (Drosophila) 

White 

1377917_at 

— 

— 

Green 

1377962_at 

Mybl2 

myeloblastosis  oneogene-like  2 

White 

1377967_at 

Cdtl 

ehromatin  lieensing  and  DNA  replieation 
faetor  1 

White 

1378056_at 

Gmnn 

geminin 

White 

1378334_a_at 

— 

— 

Yellow 

1378507_at 

— 

— 

White 

1378537_at 

— 

— 

White 

1378600_at 

Tafld 

TATA  box  binding  protein  (TBP)-assoeiated 
faetor,  RNA  polymerase  I,  D,  41kDa 

White 

1378640_at 

Uhrfl 

ubiquitin-like  with  PHD  and  ring  linger 
domains  1 

White 

1378766_at 

— 

— 

Yellow 

1378906_at 

RGD  1559690 

similar  to  hypothetieal  protein  PEJ25416 

White 

137903 l_at 

Gca 

granealcin 

Yellow 

1379073_at 

Slc25a23 

solute  earrier  family  25  (mitoehondrial 
earrier;  phosphate  carrier),  member  23 

Yellow 

1379076_at 

— 

— 

Yellow 
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BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1379393_at 

Vill 

villin  1 

Yellow 

1379607_at 

— 

— 

Yellow 

1379608_at 

RGD 15600 10 

RGD  15600 10 

White 

1379673_at 

Uaplll 

UDP-N-acteylglucosamine 
pyrophosphorylase  1-like  1 

Yellow 

1379696_at 

— 

— 

Yellow 

1379778_at 

RGD  1307983 

similar  to  HSPC043  protein 

White 

1379815_at 

LOC683733 

similar  to  Transeription  factor  7-like  2 
(HMG  box  transcription  factor  4)  (T-cel 

Yellow 

1379904_at 

Trappc6a 

trafficking  protein  particle  complex  6A 

Yellow 

1379957_at 

Slfn8 

schlafen  8 

White 

1379967_at 

Zfp367 

zinc  finger  protein  367 

White 

1379980_at 

— 

— 

Yellow 

1380013_at 

Pnpla3 

patatin-like  phospholipase  domain 
containing  3 

Yellow 

1380182_at 

RGD1563437 

similar  to  KIAA1217 

Yellow 

1380243_at 

RGD  1304693 

similar  to  CG14803-PA 

White 

1380277_at 

RadSlapl 

RAD5 1  associated  protein  1 

White 

1380320_at 

Lin54 

lin-54  homolog  (C.  elegans) 

White 

1380598_at 

— 

— 

Yellow 

1380883_at 

— 

— 

White 

1381006_at 

Hgfac 

hepatocyte  growth  factor  activator 

Yellow 

1381130_at 

Mcm8 

minichromosome  maintenance  complex 
component  8 

White 

1381574_at 

Tmeml95 

transmembrane  protein  195 

Yellow 

1381871_at 

— 

— 

Yellow 

1381982_at 

Uaplll 

UDP-N-acteylglucosamine 
pyrophosphorylase  1-like  1 

Yellow 

1381983_at 

Denndla 

DENN/MADD  domain  containing  lA 

Yellow 

1382060_at 

— 

— 

Yellow 

1382218_at 

RGD1305807 

hypothetical  LOC298077 

Yellow 

13823  ll_at 

Tifa 

TRAF-interacting  protein  with  forkhead- 
associated  domain 

Yellow 

138243 l_at 

— 

— 

Yellow 

1382659_at 

Pla2rl 

phospholipase  A2  receptor  1 

Yellow 

1382830_at 

Suv391i2 

suppressor  of  variegation  3-9  homolog  2 
(Drosophila) 

White 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1383103_at 

Faml72a 

family  with  sequence  similarity  172, 
member  A 

Yellow 

1383188_at 

— 

— 

Yellow 

1383212_at 

— 

— 

Yellow 

1383287_at 

Fubpl 

far  upstream  element  (FUSE)  binding 
protein  1 

Green 

1383458_at 

— 

— 

Yellow 

1383477_at 

Uchl5 

ubiquitin  carboxyl-terminal  hydrolase  L5 

Green 

1383497_at 

— 

— 

Yellow 

1383574_at 

Manlal 

mannosidase,  alpha,  class  lA,  member  1 

Yellow 

1383604_at 

— 

— 

Yellow 

1383637_at 

— 

— 

Yellow 

1383640_at 

— 

— 

Yellow 

1383684_at 

As  fib 

ASFl  anti-sileneing  function  1  homolog  B 
(S.  cerevisiae) 

White 

1383721_at 

Fzd8 

frizzled  homolog  8  (Drosophila) 

Green 

1383728_at 

LOC686240 

similar  to  NMD  A  reeeptor  regulated  1-like 

Green 

1383770_at 

— 

— 

Yellow 

1383853_at 

Dyrk3 

dual-specifieity  tyrosine-(Y)- 
phosphorylation  regulated  kinase  3 

White 

1383860_at 

Fosl2 

Fos-like  antigen  2 

White 

1383900_at 

Qrich2 

glutamine  rich  2 

Yellow 

1383921_at 

— 

— 

White 

138393 l_at 

— 

— 

Yellow 

1383956_at 

RGD 1565709 

similar  to  ovostatin-2 

Yellow 

1384041_at 

— 

— 

White 

1384068_at 

Ckap2 

cytoskeleton  associated  protein  2 

White 

1384103_at 

RGD1561416 

similar  to  novel  protein  (HT036) 

Yellow 

1384145_at 

RGD13 10495 

similar  to  KIAA1919  protein 

Yellow 

138423  l_at 

LOC687121 

similar  to  She  SH2-domain  binding  protein  1 

White 

1384499_at 

Bucsl 

butyryl  Coenzyme  A  synthetase  1 

Yellow 

1384701_at 

— 

— 

White 

13848 18_at 

Mylk 

myosin  light  ehain  kinase 

Yellow 

1384988_at 

Fbxo5 

F-box  protein  5 

White 

1385024_at 

Cute 

eutC  copper  transporter  homolog  (E.coli) 

Green 

1385109_at 

Tafld 

TATA  box  binding  protein  (TBP)-assoeiated 

White 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1385133_at 

LOC691729 

similar  to  spermatogenesis  associated  5-like 

1 

White 

1385381_at 

— 

Yellow 

1385522_at 

Ore  11 

origin  recognition  complex,  subunit  1-like 
(yeast) 

White 

1385733_at 

Exol 

exonuelease  1 

White 

1385756_at 

Itihl 

inter-alpha  trypsin  inhibitor,  heavy  chain  1 

Yellow 

1385781_at 

Ercc61 

excision  repair  eross-complementing  rodent 
repair  deficiency  complementation  gro 

White 

1386466_at 

— 

— 

White 

1387034_at 

Pah 

phenylalanine  hydroxylase 

Yellow 

1387045_at 

Atp6v0al 

ATPase,  H-H  transporting,  lysosomal  VO 
subunit  A1 

Yellow 

1387050_s_at 

Kngl  /// 
Knglll  /// 
Kng2 

kininogen  1  III  kininogen  1-like  1  III 
kininogen  2 

Yellow 

1387307_at 

Hal 

histidine  ammonia  lyase 

Yellow 

138741 l_at 

Ptprk 

protein  tyrosine  phosphatase,  receptor  type, 

K,  extraeellular  region 

Yellow 

1387426_at 

Slc25a21 

solute  earner  family  25  (mitoehondrial 
oxodiearboxylate  carrier),  member  21 

Yellow 

1387457_at 

Duspl2 

dual  specifieity  phosphatase  12 

Green 

1387530  a  at 

Posh  ///  Eosl2 

EBJ  osteosarcoma  oncogene  B  III  fos-like 
antigen  2 

White 

1387599_a_at 

Nqol 

NAD(P)H  dehydrogenase,  quinone  1 

White 

1387626_at 

Dck 

deoxycytidine  kinase 

White 

1387673  a  at 

Anxa6 

annexin  A6 

Yellow 

1387725_at 

Gulo 

gulonolactone  (E-)  oxidase 

Yellow 

1387795_at 

Pola2 

polymerase  (DNA  direeted),  alpha  2 

White 

1387872_at 

Hnrnpal 

heterogeneous  nuclear  ribonucleoprotein  A1 

Green 

1387907_at 

Itprl 

inositol  1,4,5-triphosphate  receptor,  type  1 

Yellow 

1387924_at 

Ngef 

neuronal  guanine  nueleotide  exchange  factor 

Yellow 

1387963  a  at 

Uox 

urate  oxidase 

Yellow 

1387972_at 

Mucdhl 

muein  and  cadherin  like 

Yellow 

1388135_at 

Rpa2 

replication  protein  A2 

White 

1388143_at 

Coll8al 

eollagen,  type  XVIII,  alpha  1 

Yellow 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1388639_at 

Bcas3 

Breast  careinoma  amplified  sequenee  3 

Yellow 

1388744_at 

Mcm7 

miniehromosome  maintenanee  defieient  7 
(S.  eerevisiae) 

White 

1388825_at 

— 

— 

Yellow 

1388870_at 

Msi2 

Musashi  homolog  2  (Drosophila) 

Y  ellow 

1389102_at 

Dockl 

dedieator  of  eyto-kinesis  1 

Yellow 

1389136_at 

Prkrir 

protein-kinase,  interferon-indueible  double 
stranded  RNA  dependent  inhibitor,  re 

White 

1389150_at 

— 

— 

Yellow 

1389159_at 

— 

— 

Yellow 

1389248_at 

Galkl 

galaetokinase  1 

Yellow 

1389371_at 

— 

— 

Yellow 

1389408_at 

Rrm2 

ribonucleotide  reductase  M2 

White 

1389421_at 

Pole 

polymerase  (DNA  directed),  epsilon 

White 

1389434_at 

Bbs9 

Bardet-Biedl  syndrome  9 

Yellow 

1389439_at 

— 

— 

Green 

1389440_at 

— 

— 

White 

1389448_at 

4-Sep 

septin  4 

Y  ellow 

1389457_at 

Mybl2 

myeloblastosis  oncogene-like  2 

White 

13895  lOat 

Lyar 

Lyl  antibody  reactive  homolog  (mouse) 

White 

1389649_at 

— 

— 

Yellow 

1389725_at 

Tm7sf2 

transmembrane  7  superfamily  member  2 

Yellow 

1389852_at 

Cenpm 

centromere  protein  M 

White 

1390389_at 

— 

— 

Yellow 

1390404_at 

Lama2 

laminin,  alpha  2 

Yellow 

1390432_at 

Tmco3 

transmembrane  and  coiled-coil  domains  3 

Yellow 

1 39048  l_a_at 

Ube2t 

ubiquitin-conjugating  enzyme  E2T 
(putative) 

White 

1390557_at 

Gca 

grancalcin 

Yellow 

1390623_at 

Gins3 

GINS  complex  subunit  3  (Psf3  homolog) 

White 

1390650_at 

Nup85 

nucleoporin  85kDa 

White 

1390729_x_at 

RGD 1305441 

similar  to  RIKEN  cDNA  5730453116 

Green 

1390823_at 

Skp2 

S-phase  kinase-associated  protein  2  (p45) 

White 

1390856_at 

— 

— 

Yellow 

1390982_at 

Chaflb 

chromatin  assembly  factor  1 ,  subunit  B 

(p60) 

White 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1391526_at 

Kifl6b 

kinesin  family  member  16B 

Y  ellow 

1391538_at 

Cideb 

eell  death-indueing  DNA  fragmentation 
faetor,  alpha  subunit-like  effeetor  B 

Yellow 

1391544_at 

— 

— 

Yellow 

139161  l_at 

Kifl2 

kinesin  family  member  12 

Yellow 

1391621_at 

— 

— 

Yellow 

1391871_at 

Stardl3 

StAR-related  lipid  transfer  (START)  domain 
eontaining  13 

Yellow 

1391890_at 

Bcas3 

breast  eareinoma  amplified  sequenee  3 

Yellow 

1391897_at 

Rabgapll 

RAB  GTPase  aetivating  protein  1-like 

Yellow 

1391916_at 

— 

— 

Yellow 

1392044_at 

— 

— 

Yellow 

1392057_at 

LOC500893 

similar  to  GLI-Kruppel  family  member 

GLI4 

Blue 

1 39206  l_at 

Mem  10 

miniehromosome  maintenanee  eomplex 
eomponent  10 

White 

1392064_at 

Dlxl 

distal-less  homeobox  1 

White 

1392128_at 

Atad5 

ATPase  family,  AAA  domain  eontaining  5 

White 

1392440_at 

— 

— 

White 

1392554_a_at 

— 

— 

White 

1392653_at 

— 

— 

Yellow 

1392720_at 

Cyp4fl7 

eytoehrome  P450,  family  4,  subfamily  f, 
polypeptide  17 

Yellow 

1392800_at 

RGD 1308 165 

similar  to  hypothetieal  protein  MGC 17337 

Blue 

1392825_at 

RGD 1559600 

RGD  1559600 

Yellow 

1392852_at 

Smocl 

SPARC  related  modular  ealeium  binding  1 

Yellow 

1392870_at 

— 

— 

Yellow 

1393036_at 

Psrcl 

proline/serine -rieh  eoiled-eoil  1 

White 

1393123_at 

C8g 

eomplement  eomponent  8,  gamma 
polypeptide 

Yellow 

1393228_at 

Hpdl 

4-hydroxyphenylpyruvate  dioxygenase-like 

White 

1393230_s_at 

RGD  1308772 

Similar  to  KIAA0564  protein 

Yellow 

1393260_at 

LOC360228 

WDNMl  homolog 

Yellow 

1393373_at 

Opn3 

opsin  3 

White 

1393561_at 

— 

— 

Yellow 

1393803_at 

Casp8ap2 

easpase  8  assoeiated  protein  2 

White 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1393982_at 

Pole2 

polymerase  (DNA  directed),  epsilon  2  (p59 
subunit) 

White 

1394086_at 

Senp7 

SUMOl/sentrin  specific  protease  7 

Yellow 

1394371_at 

Nufipl 

nuclear  fragile  X  mental  retardation  protein 
interacting  protein  1 

Blue 

1394419_at 

Arhgapl la 

Rho  GTPase  activating  protein  1 1 A 

White 

1394457_at 

— 

— 

White 

1395132_at 

Utrn 

Utrophin 

Y  ellow 

1395203_at 

— 

— 

Yellow 

1395274_at 

Dst 

dystonin 

Yellow 

1395347_at 

Arsa 

arylsulfatase  A 

Yellow 

1395644_at 

Abca8 

ATP -binding  cassette,  sub-family  A 
(ABCl),  member  8 

Yellow 

1395944_at 

RGD13 10778 

similar  to  Putative  protein  C21orf45 

White 

1395957_at 

— 

— 

Yellow 

1396236_at 

— 

— 

White 

1396323_at 

LOC681383 

similar  to  Protein  ClOorfl  1  homolog 

Yellow 

1397286_at 

— 

— 

Yellow 

1397406_at 

— 

— 

Yellow 

1397409_s_at 

Weel 

wee  1  homolog  (S.  pombe) 

White 

1397556_at 

Natl  3 

N-acetyltransferase  13 

Green 

1397564_at 

— 

— 

White 

1397698_at 

— 

— 

Yellow 

1397704_at 

LOC691979 

similar  to  N-acetyltransferase  ESC02 
(Establishment  of  cohesion  1  homolog  2) 

White 

1397915_at 

— 

— 

White 

1398397_at 

— 

— 

Yellow 

1398439_a_at 

Orc61 

origin  recognition  complex,  subunit  6  like 
(yeast) 

White 

1398546_at 

— 

— 

White 

1398598_at 

Dst 

dystonin 

Yellow 

1398706_at 

— 

— 

White 

137413 l_at 

— 

— 

Yellow 

1369630_at 

Adk 

adenosine  kinase 

Yellow 

1371862_at 

Rrml 

ribonucleotide  reductase  Ml 

White 

1376485_at 

Traip 

TRAE-interacting  protein 

White 
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APPENDIX  3  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM 


Affymetix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1383063  a  at 

Yellow 

1385132_at 

— 

— 

White 

1388526_at 

Gstzl 

glutathione  transferase  zeta  1 

Yellow 

1393848_at 

Rrm2 

ribonueleotide  reduetase  M2 

White 

1374087_at 

Slc48al 

solute  earrier  family  48  (heme  transporter), 
member  1 

White 

1379277_at 

— 

— 

Yellow 

1383002_at 

PhfZOll 

PHD  linger  protein  20-like  1 

Yellow 

1387082_at 

Fetub 

fetuin  B 

Yellow 

1387221_at 

Gchl 

GTP  eyelohydrolase  1 

Green 

1391077_at 

Clspn 

elaspin  homolog  (Xenopus  laevis) 

White 

1393844_at 

Mlflip 

myeloid  leukemia  faetor  1  interaeting 
protein 

White 

1379386_at 

Gtsel 

G-2  and  S-phase  expressed  1 

White 
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APPENDIX  4:  PROBE  SETS  DIFFERENTIALLY  EXPRESSED  BY  CADMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

136763  l_at 

Ctgf 

connective  tissue  growth  factor 

Green 

1367650_at 

Tinagll 

tubulointerstitial  nephritis  antigen-like  1 

Green 

1367655_at 

TmsblO 

thymosin,  beta  10 

Green 

1367661  at 

S100a6 

SI 00  calcium  binding  protein  A6 

Green 

1367729_at 

Oat 

ornithine  aminotransferase  (gyrate  atrophy) 

Yellow 

1367764_at 

Ccngl 

cyclin  G1 

Green 

1367784_a_at 

Clu 

clusterin 

Green 

1367800_at 

Plat 

plasminogen  activator,  tissue 

Green 

1367974_at 

Anxa3 

annexin  A3 

Green 

1368072_at 

Btg3 

B-cell  translocation  gene  3 

Green 

1368162_at 

Cst6 

cystatin  E/M 

Green 

1368167_at 

Ctse 

cathepsin  E 

Green 

1368205_at 

Cfi 

complement  factor  I 

Yellow 

1368247_at 

Hspala/// 

Hspalb 

heat  shock  70kD  protein  1 A  III  heat  shock 
70kD  protein  IB  (mapped) 

Green 

1368269_at 

Lgals4 

lectin,  galactoside-binding,  soluble,  4 

Green 

1368290_at 

Cyr61 

cysteine-rich,  angiogenic  inducer,  61 

Green 

1368332_at 

Gbp2 

guanylate  binding  protein  2 

Green 

1368335_at 

Apoal 

apolipoprotein  A-I 

Green 

1368339_at 

SlOOg 

SI 00  calcium  binding  protein  G 

Green 

1368379_at 

Scarb2 

scavenger  receptor  class  B,  member  2 

Green 

1368380_at 

Vtn 

vitronectin 

Yellow 

1368490_at 

Cdl4 

CD  14  molecule 

Green 

1368538_at 

Exoc7 

exocyst  complex  component  7 

Green 

1368905_at 

Ces2 

carboxylesterase  2  (intestine,  liver) 

Green 

1368947_at 

Gadd45a 

growth  arrest  and  DNA-damage-inducible, 
alpha 

Green 

1369202_at 

Mx2 

myxovirus  (influenza  virus)  resistance  2 

Green 

1369286_at 

Proc 

protein  C 

Yellow 

1369590_a_at 

Ddit3 

DNA-damage  inducible  transcript  3 

Green 

1369701_at 

Lipc 

lipase,  hepatic 

Yellow 

1370057_at 

Csrpl 

cysteine  and  glycine -rich  protein  1 

Green 

1370073_at 

Dnajc3 

DnaJ  (Hsp40)  homolog,  subfamily  C, 
member  3 

Green 

1370153_at 

Gdfl5 

growth  differentiation  factor  15 

Green 

1370174_at 

Ppplrl5a 

protein  phosphatase  1 ,  regulatory  (inhibitor) 
subunit  15A 

Green 

1370223_at 

Arfrpl 

ADP-ribosylation  factor  related  protein  1 

Green 
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APPENDIX  4  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CADMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1370259_a_at 

Pthrl 

parathyroid  hormone  reeeptor  1 

Green 

13703 13_at 

Acot7 

acyl-CoA  thioesterase  7 

Green 

1370348_at 

Ninjl 

ninjurin  1 

Green 

1370361_at 

Cgrefl 

cell  growth  regulator  with  EF  hand  domain  1 

Green 

1370387_at 

Cyp3a9 

cytochrome  P450,  family  3,  subfamily  a, 
polypeptide  9 

Green 

1370399_at 

Cyp4bl 

cytochrome  P450,  family  4,  subfamily  b, 
polypeptide  1 

Green 

1 37040  l_at 

Ly6b 

lymphocyte  antigen  6  complex,  locus  B 

Green 

1370464_at 

Abcbl 

ATP-binding  cassette,  sub-family  B 
(MDR/TAP),  member  1 

Green 

1370583_s_at 

Abcbl  /// 
Abcblb 

ATP-binding  cassette,  sub-family  B 
(MDR/TAP),  member  1  III  ATP-binding 
cassette. 

Green 

1370694_at 

Trib3 

tribbles  homolog  3  (Drosophila) 

Green 

1370803_at 

Zwint 

ZWIO  interactor 

Green 

1370912_at 

Hspalb 

heat  shock  70kD  protein  IB  (mapped) 

Green 

1370988_at 

— 

— 

Green 

1371015_at 

Mxl 

myxovirus  (influenza  virus)  resistance  1 

Green 

137121 l_a_at 

Nrgl 

neuregulin  1 

Green 

1371212_at 

Nrgl 

neuregulin  1 

Green 

1371455_at 

Pmml 

phosphomannomutase  1 

Green 

1371689_at 

— 

— 

Yellow 

1371692_at 

Mlltl  1 

myeloid/lymphoid  or  mixed-lineage 
leukemia  (trithorax  homolog.  Drosophila) 

Green 

1371703_at 

Ahnak 

AHNAK  nucleoprotein 

Green 

1371735_at 

— 

— 

Green 

1371774_at 

Sail 

spermidine/spermine  N1 -acetyl  transferase  1 

Green 

1371785_at 

Tnfrsfl2a 

tumor  necrosis  factor  receptor  superfamily, 
member  12a 

Green 

1372216_at 

Cyb561dl 

cytochrome  b-5  61  domain  containing  1 

Green 

1372240_at 

Sgca 

sarcoglycan,  alpha  (dystrophin-associated 
glycoprotein) 

Green 

1372306_at 

Ethel 

ethylmalonic  encephalopathy  1 

Green 

1372476_at 

Fads3 

fatty  acid  desaturase  3 

Green 

1372585_at 

RGD 1566254 

RGD  1566254 

Green 

1372604_at 

RGD 1309808 

similar  to  apolipoprotein  E2;  apolipoprotein 

Green 

L-II 
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APPENDIX  4  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CADMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1372762_at 

Hddc3 

HD  domain  containing  3 

Green 

1388190_at 

Apob 

apolipoprotein  B  (including  Ag(x)  antigen) 

Yellow 

1388198_at 

Nupll 

nucleoporin  like  1 

White 

1372809_at 

LOC290595 

hypothetieal  gene  supported  by  AF 152002 

Green 

1372870_at 

Kdelr3 

KDEL  (Eys-Asp-Glu-Eeu)  endoplasmie 
retieulum  protein  retention  reeeptor  3 

Green 

1373490_at 

Gmfg 

glia  maturation  faetor,  gamma 

Green 

1373504_at 

Gliprl 

GEl  pathogenesis-related  1 

Green 

1373606_at 

Zfand2b 

zine  finger,  AN  1  type  domain  2B 

Green 

1374146_at 

Mad212 

MAD2  mitotie  arrest  defieient-like  2  (yeast) 

Green 

1374345_at 

— 

___ 

Green 

1374610_at 

— 

— 

Green 

1375170_at 

SlOOall 

SI 00  ealeium  binding  protein  A1 1 
(ealizzarin) 

Green 

1375346_at 

RGD 1563941 

similar  to  hypothetieal  protein  FEJ20010 

Green 

1375684_at 

Neul 

sialidase  1  (lysosomal  sialidase) 

Green 

1375906_at 

LOC683034 

hypothetieal  protein  LOC683034 

Green 

1375920_at 

Plekhol 

pleekstrin  homology  domain  eontaining, 
family  0  member  1 

Green 

13763 19_at 

Sema3c 

sema  domain,  immunoglobulin  domain  (Ig), 
short  basie  domain,  seereted 

Green 

1376500_at 

Fbxo23 

F-box  only  protein  23 

Green 

1377336_at 

Sema3b 

sema  domain,  immunoglobulin  domain  (Ig), 
short  basie  domain,  seereted,  (semaphor 

Green 

1377456_at 

— 

— 

Yellow 

1377739_at 

Gngl2 

guanine  nueleotide  binding  protein  (G 
protein),  gamma  12 

Green 

1377752_at 

Emp2 

epithelial  membrane  protein  2 

Green 

1378387_at 

Tinf2 

TERFl  (TRFl)-interaeting  nuelear  faetor  2 

Green 

1378570_at 

— 

— 

Yellow 

1378883_at 

— 

— 

Yellow 

1378896_at 

Slc30a2 

solute  earner  family  30  (zine  transporter), 
member  2 

Green 

1378910_at 

— 

— 

Green 

1379318_at 

LOC498662 

similar  to  RIKEN  eDNA  2610019F03 

Green 

1379420_at 

RGD  1565002 

similar  to  Dehydrogenase/reduetase  SDR 
family  member  7  preeursor  (Retinal  short- 

Yellow 

1379805_at 

Slc41a2 

solute  earrier  family  41,  member  2 

Green 
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APPENDIX  4  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CADMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1380470_x_at 

— 

— 

Green 

138061 l_at 

Fkbp5 

FK506  binding  protein  5 

Green 

1381082_at 

RGD 1564357 

RGD1564357 

Green 

1381428_a_at 

— 

— 

Yellow 

1381453_at 

— 

— 

Yellow 

1381940_at 

— 

___ 

Green 

1382028_at 

Pexl9 

peroxisomal  biogenesis  faetor  19 

Yellow 

1382137_at 

Abhd3 

abhydrolase  domain  eontaining  3 

Green 

1382475_at 

— 

... 

Green 

1382554_at 

C8a 

complement  component  8,  alpha  polypeptide 

Yellow 

1382882_x_at 

— 

... 

Green 

1383248_at 

Fmo5 

flavin  containing  monooxygenase  5 

Green 

1383302_at 

Dnajbl 

DnaJ  (Hsp40)  homolog,  subfamily  B, 
member  1 

Green 

1383369_at 

Trim26 

tripartite  motif-containing  26 

Green 

1383485_at 

Mdm2 

Mdm2  p53  binding  protein  homolog 
(mouse) 

Green 

1383666_at 

Ptrhl 

peptidyl-tRNA  hydrolase  1  homolog  (S. 
cerevisiae) 

Green 

1383702_at 

— 

— 

Yellow 

1383797_a_at 

Ptrhl 

peptidyl-tRNA  hydrolase  1  homolog  (S. 
cerevisiae) 

Green 

1383874_at 

RGD1560812 

RGD1560812 

Green 

1383999_at 

— 

... 

Yellow 

1384602_at 

Gprl07 

G  protein-coupled  receptor  107 

Green 

1385001_at 

Gsdmd 

gasdermin  D 

Green 

1385089_at 

Grhll 

grainyhead-like  1  (Drosophila) 

Green 

13851 17_at 

— 

... 

Green 

1385252_at 

Trim6 

tripartite  motif-containing  6 

Green 

1385269_s_at 

Ccdc93 

coiled-coil  domain  containing  93 

Green 

1385724_at 

LOC682874 

hypothetical  protein  LOC682874 

Green 

1385801_at 

Dnajcl8 

DnaJ  (Hsp40)  homolog,  subfamily  C, 
member  18 

Green 

1385899_at 

Triml6 

tripartite  motif-containing  16 

Green 

1386186_s_at 

— 

... 

Green 

1386879_at 

Lgals3 

lectin,  galactoside-binding,  soluble,  3 

Green 

1386890_at 

SlOOalO 

SI 00  calcium  binding  protein  AlO 

Green 

1386901_at 

Cd36 

CD36  molecule  (thrombospondin  receptor) 

Y  ellow 
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APPENDIX  4  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CADMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1387291_at 

Itih3 

inter-alpha  trypsin  inhibitor,  heavy  ehain  3 

Yellow 

13873 16_at 

Cxcll 

ehemokine  (C-X-C  motif)  ligand  1 
(melanoma  growth  stimulating  activity, 
alpha) 

Green 

1387391_at 

Cdknla 

cyclin-dependent  kinase  inhibitor  lA  (p21, 
Cipl) 

Green 

138745 l_at 

Dcbld2 

discoidin,  CUB  and  LCCL  domain 
containing  2 

Green 

1387663_at 

Gmfb 

glia  maturation  factor,  beta 

Green 

13878 18_at 

Casp4 

caspase  4,  apoptosis-related  cysteine 
peptidase 

Green 

1387856_at 

Cnn3 

calponin  3,  acidic 

Green 

1388130_at 

Zyx 

zyxin 

Green 

1388146_at 

— 

— 

Yellow 

1388453_at 

Myadm 

myeloid-associated  differentiation  marker 

Green 

1388460_at 

Capg 

capping  protein  (actin  fdament),  gelsolin¬ 
like 

Green 

1388574_at 

Wars 

tryptophanyl-tRNA  synthetase 

Green 

1388663_at 

LOC687118 

similar  to  death  effector  domain-containing 
DNA  binding  protein  2 

Green 

1388722_at 

Dnajbl 

DnaJ  (Hsp40)  homolog,  subfamily  B, 
member  1 

Green 

1388763_at 

LOC690976 

similar  to  Calponin-2  (Calponin  H2,  smooth 
muscle)  (Neutral  calponin) 

Green 

1388792_at 

Gadd45g 

growth  arrest  and  DNA-damage-inducible, 
gamma 

Green 

1388924_at 

Angptl4 

angiopoietin-like  4 

Green 

1389068_at 

Wipi2 

WD  repeat  domain,  phosphoinositide 
interacting  2 

Green 

13891  lOat 

— 

... 

Green 

1389166_at 

Cib2 

calcium  and  integrin  binding  family  member 

2 

Green 

1389402_at 

Axudl 

AXINl  up-regulated  1 

Green 

1389650_at 

— 

... 

Yellow 

1389732_at 

Dram 

damage-regulated  autophagy  modulator 

Green 

1389990_at 

— 

... 

Green 

1390026_at 

Bag3 

Bcl2-associated  athanogene  3 

Green 

1390042_at 

Tmeml40 

transmembrane  protein  140 

Yellow 

1390391  at 

— 

— 

Green 
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APPENDIX  4  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CADMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1 39069  l_at 

RGD1561551 

Similar  to  Hypothetical  protein  MGC75664 

Green 

1390770_at 

— 

— 

Green 

1391028_at 

— 

___ 

Green 

1391571_at 

— 

___ 

Yellow 

1392644_s_at 

LOC498662 

similar  to  RIKEN  cDNA  2610019F03 

Green 

1393025_at 

— 

___ 

Green 

1393134_at 

Sephsl 

selenophosphate  synthetase  1 

Yellow 

1393584_at 

TnfrsfZl 

tumor  necrosis  factor  receptor  superfamily, 
member  21 

Green 

1394363_at 

Riok3 

RIO  kinase  3  (yeast) 

Green 

1394923_at 

LOC689619 

similar  to  nuclear  receptor  interacting 
protein  2 

Green 

1395329_at 

Gmppb 

GDP-mannose  pyrophosphorylase  B 

Green 

1395366_at 

— 

— 

Green 

1395572_at 

Dnaja4 

DnaJ  (Hsp40)  homolog,  subfamily  A, 
member  4 

Green 

1395586_at 

— 

— 

Yellow 

1395699_at 

Riok3 

RIO  kinase  3  (yeast) 

Green 

1396228_at 

RGD1308958 

Similar  to  chromosome  9  open  reading 
frame  5 

Green 

1396383_at 

Aldh3bl 

aldehyde  dehydrogenase  3  family,  member 

B1 

Green 

1 39656  l_x_at 

Piga 

phosphatidylinositol  glycan  anchor 
biosynthesis,  class  A 

Green 

1397372_at 

Mobkllb 

MOBl,  Mps  One  Binder  kinase  activator¬ 
like  IB  (yeast) 

Green 

1398245_at 

Sncg 

synuclein,  gamma  (breast  cancer-specific 
protein  1) 

Green 

1398264_at 

Slc30a2 

solute  carrier  family  30  (zinc  transporter), 
member  2 

Green 

1398343_at 

LOC498996 

similar  to  DnaJ  (Hsp40)  homolog,  subfamily 
A,  member  4 

Green 

1398376_at 

— 

... 

Green 

1398600_at 

— 

... 

Green 

1398771_at 

Slc3a2 

solute  carrier  family  3  (activators  of  dibasic 
and  neutral  amino  acid  transport) 

Green 

1370244_at 

Ctsll 

cathepsin  LI 

Green 

1380668_at 

MGC72974 

Hypothetical  LOC3 16976 

Green 

1387282_at 

Hspb8 

heat  shock  protein  8 

Green 
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APPENDIX  4  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CADMIUM 


1390662 

at 

Trim24 

tripartite  motif-containing  24 

Green 

1372729_ 

at 

Procr 

protein  C  receptor,  endothelial 

Green 

1379854 

at 

Abhd5 

abhydrolase  domain  containing  5 

Green 

138701 1_ 

at 

Lcn2 

lipocalin  2 

Green 

1390455 

at 

Abhd2 

abhydrolase  domain  containing  2 

Green 
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APPENDIX  5:  PROBE  SETS  DIFFERENTIALLY  EXPRESSED  BY  NICKEL, 

CHROMIUM,  AND  CADMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1367850_at 

Fcgr2a  — 
LOC498276  - 
-  LOC498277 

Fc  fragment  of  IgG,  low  affinity  Ila,  receptor 
(CD32)  —  Fc  gamma  receptor  II  b 

Green 

1368342_at 

Ampd3 

adenosine  monophosphate  deaminase  3 

Green 

1368363_at 

Klf5 

Kruppel-like  factor  5 

Green 

1368550_at 

Foxql 

forkhead  box  Q1 

Green 

1369050_at 

Pik3c2g 

phosphoinositide-3-kinase,  class  2,  gamma 
polypeptide 

Yellow 

13691  Hat 

Fabpl 

fatty  acid  binding  protein  1,  liver 

Yellow 

1369268_at 

Atfi 

activating  transcription  factor  3 

Green 

1370080_at 

Hmoxl 

heme  oxygenase  (decycling)  1 

Green 

1370491_a_at 

Hdc 

histidine  decarboxylase 

Yellow 

1370725_a_at 

G6pc 

glucose-6-phosphatase,  catalytic  subunit 

Yellow 

1371109_at 

C8b 

complement  component  8,  beta  polypeptide 

Yellow 

1371237_a_at 

Mtla 

metallothionein  la 

Green 

1371527_at 

Empl 

epithelial  membrane  protein  1 

Green 

1372727_at 

— 

— 

Green 

1373767_at 

Zfand2a 

zinc  linger,  AN  1  -type  domain  2A 

Green 

1374886_at 

Bcsll 

BCSl-like  (yeast) 

Blue 

1380229_at 

Maff 

v-maf  musculoaponeurotic  fibrosarcoma 
oncogene  homolog  F  (avian) 

Blue 

1380699_at 

— 

— 

Green 

1383253_at 

— 

— 

Green 

1383425_at 

C5 

complement  component  5 

Yellow 

1383632_at 

— 

— 

Green 

1385961_at 

Klf5 

Kruppel-like  factor  5 

Green 

1387228_at 

Slc2a2 

solute  carrier  family  2  (facilitated  glucose 
transporter),  member  2 

Yellow 

1388153_at 

Acsll 

acyl-CoA  synthetase  long-chain  family 
member  1 

Yellow 

1388271_at 

Mt2A 

metallothionein  2A 

Green 

1 38938  l_at 

Sqstml 

sequestosome  1 

Green 

1391806_at 

LOC498793 

similar  to  inter-alpha-inhibitor  H2  chain 

Yellow 

1394039_at 

Klf5 

Kruppel-like  factor  5 

Green 

13941 12_at 

Haol 

hydroxyacid  oxidase  1 

Yellow 
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APPENDIX  5  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 
BY  NICKEL,  CHROMIUM,  AND  CADMIUM 


Affymetrix 
ID _ 

1398260_a_at 

1398246  s  at 


Gene  Symbol  Gene  Title 


Serpindl 

Fcgr2a 

LOC498276 


serine  (or  cysteine)  peptidase  inhibitor,  clade 
D,  member  1 

Fc  fragment  of  IgG,  low  affinity  Ila,  receptor 
(CD32)  —  Fc  gamma  receptor  II  b 


Cluster 

Yellow 

Green 
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APPENDIX  6:  PROBE  SETS  DIFFERENTIALLY  EXPRESSED  BY  NICKEL 

AND  CADMIUM 


Affymetrix 

ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1367584_at 

Anxa2 

annexin  A2 

Green 

1367593_at 

Sepwl 

selenoprotein  W,  1 

Green 

1367795_at 

Ifrdl 

interferon-related  developmental  regulator  1 

Blue 

1368146_at 

Duspl 

dual  speeifieity  phosphatase  1 

Green 

1368147_at 

Duspl 

dual  speeifieity  phosphatase  1 

Green 

1368940_at 

P2ry2 

purinergie  reeeptor  P2Y,  G-protein  eoupled  2 

Green 

1369099_at 

Slc30al 

solute  earner  family  30  (zinc  transporter), 
member  1 

Green 

1369736_at 

Empl 

epithelial  membrane  protein  1 

Green 

1370000_at 

Nucb2 

nucleobindin  2 

Blue 

1370026_at 

Cryab 

crystallin,  alpha  B 

Green 

1370177_at 

PVR 

poliovirus  receptor 

Green 

1371360_at 

Ndrgl 

N-myc  downstream  regulated  gene  1 

Blue 

1372103_at 

LOC498996 

similar  to  DnaJ  (Hsp40)  homolog,  subfamily 
A,  member  4 

Green 

1372293_at 

RGD 1562059 

Similar  to  RIKEN  cDNA  1 1 10038E21 

Blue 

1373677_at 

Slc39al0 

solute  carrier  family  39  (zinc  transporter), 
member  10 

Yellow 

1374366_at 

Slc39a4 

solute  carrier  family  39  (zinc  transporter), 
member  4 

Yellow 

1374404_at 

Jun 

Jun  oncogene 

Green 

1376513  a  at 

— 

— 

Green 

1377626_at 

LOC690768 

hypothetical  protein  EOC690768 

Yellow 

1379934_at 

— 

— 

Blue 

1381335_at 

— 

— 

Green 

1381533_at 

Rndl 

Rho  family  GTPase  1 

Green 

1383193_at 

— 

... 

Green 

1383744_at 

— 

... 

Green 

1385464_at 

Foxql 

Eorkhead  box  Q1 

Green 

1386944_a_at 

G6pc 

glucose-6-phosphatase,  catalytic  subunit 

Yellow 

1388267_a_at 

Mtla 

metallothionein  la 

Green 

1388721_at 

Hspb8 

heat  shock  protein  8 

Green 

1388900_at 

RGD1566118 

RGD1566118 

Green 

1389528_s_at 

Jun 

Jun  oncogene 

Green 
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APPENDIX  6  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY 
EXPRESSED  BY  NICKEL  AND  CADMIUM 


Affymetrix  ID 

Gene 

Symbol 

Gene  Title 

Cluster 

1391481_at 

— 

— 

Green 

1392579_at 

Obfc2a 

oligonucleotide/oligosaccharide-binding  fold 
containing  2A 

Green 

1392627_x_at 

— 

— 

Green 

1398296_at 

Gdel 

glycerophosphodiester  phosphodiesterase  1 

Blue 

1390287_at 

Bin2  — 
LOC685015 

bridging  integrator  2  —  hypothetical  protein 
LOC685015 

Yellow 
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APPENDIX  7:  PROBE  SETS  DIFFERENTIALLY  EXPRESSED  BY  NICKEL 

AND  CHROMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1367598_at 
1367909_at 
1368180  s  at 

1368278_at 

1368627_at 

Ttr 

Dcxr 

Gsta2  /// 

LOC494499 

Lgals2 

Rgn 

transthyretin 

dicarbonyl  L-xylulose  reductase 
glutathione  S-transferase  A2  III  LOC494499 
protein 

lectin,  galactoside -binding,  soluble  2 
regucalcin  (senescence  marker  protein-30) 

Yellow 

Yellow 

Blue 

Yellow 

Yellow 

1369499_at 

Tyms 

thymidylate  synthetase 

White 

1370019_at 

Sultlal 

sulfotransferase  family,  cytosolic,  lA, 
phenol-preferring,  member  1 

Yellow 

1370264_at 

Synel 

spectrin  repeat  containing,  nuclear  envelope 

1 

3  -hydroxy-3  -methylglutaryl-Coenzyme  A 
synthase  2  (mitochondrial) 

Yellow 

13703  lOat 

Hmgcs2 

Yellow 

1370902_at 

Akrlb8 

aldo-keto  reductase  family  1,  member  B8 

Blue 

1370939_at 

Acsll 

acyl-CoA  synthetase  long-chain  family 
member  1 

Yellow 

1372006_at 

— 

— 

Blue 

1372101_at 

Ppap2b 

phosphatidic  acid  phosphatase  type  2B 

Yellow 

1372398_at 

Kctdl5 

potassium  channel  tetramerisation  domain 
containing  15 

Blue 

1373035_at 

— 

— 

Blue 

1374615_at 

— 

___ 

Yellow 

1374758_at 

— 

___ 

Yellow 

1376640_at 

Faml71al 

Family  with  sequence  similarity  171, 
member  A1 

Yellow 

1377758_at 

Hsdl7bl3 

hydroxysteroid  (17-beta)  dehydrogenase  13 

Yellow 

1377878_at 

Fgfbp3 

fibroblast  growth  factor  binding  protein  3 

Green 

1378133_at 

Slc7all 

solute  carrier  family  7  (cationic  amino  acid 
transporter,  y-H  system),  member  1 1 

Blue 

1379240_at 

— 

... 

Yellow 

1380610_at 

— 

— 

Blue 

1380682_at 

Mex3b 

mex3  homolog  B  (C.  elegans) 

Blue 

1383472_at 

Aldhlbl 

aldehyde  dehydrogenase  1  family,  member 

B1 

Yellow 

1383578_at 

Rad51 

RAD5 1  homolog  (RecA  homolog,  E.  coli) 

(S.  cerevisiae) 

White 

1383603_at 

— 

... 

Yellow 

1384505_at 

— 

... 

Yellow 

1385569  at 

— 

— 

Blue 
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APPENDIX  7  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  NICKEL  AND  CHROMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1387808_at 

Slc7a7 

solute  carrier  family  7  (cationic  amino  acid 
transporter,  y+  system),  member  7 

Yellow 

1387809_at 

Map2k6 

mitogen-activated  protein  kinase  kinase  6 

Yellow 

1388340_at 

Ns5atp9 

NS5A  (hepatitis  C  virus)  transactivated 
protein  9 

White 

1389054_at 

LOC498368 

similar  to  RIKEN  cDNA  06 10040 JO  1 

Yellow 

1389265_at 

Gbel 

glucan  (1,4-alpha-),  branching  enzyme  1 

Blue 

1389662_at 

Wnk4 

WNK  lysine  deficient  protein  kinase  4 

Yellow 

1390024_at 

— 

___ 

Yellow 

1390715_at 

Igfbpll 

insulin-like  growth  factor  binding  protein¬ 
like  1 

Yellow 

1391450_at 

Loxl2 

lysyl  oxidase-like  2 

Blue 

1392037_at 

— 

— 

Yellow 

1392384_s_at 

Akrldl 

aldo-keto  reductase  family  1,  member  D1 
(delta  4-3  -ketosteroid-5  -beta-reductase) 

Yellow 

13931 19_at 

— 

— 

Green 

1393615_at 

Depdc6 

DEP  domain  containing  6 

Yellow 

1396189_at 

— 

— 

Blue 

1396190_x_at 

— 

— 

Blue 

1398620_at 

Samd4a 

Sterile  alpha  motif  domain  containing  4 A 

Blue 

1398634_at 

LOC686032 

hypothetical  protein  EOC686032 

Yellow 

1387234_at 

Azgpl 

alpha-2-glycoprotein  1 ,  zinc -binding 

Yellow 
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APPENDIX  8:  PROBE  SETS  DIFFERENTIALLY  EXPRESSED  BY 
CHROMIUM  AND  CADMIUM 


Affymetrix 

ID 

Gene  Symbol 

Gene  Title 

Cluster 

1367647_at 

Serpinal 

serpin  peptidase  inhibitor,  elade  A  (alpha- 1 
antiproteinase,  antitrypsin),  membe 

Yellow 

1367988_at 

Cyp2c23 

eytochrome  P450,  family  2,  subfamily  c, 
polypeptide  23 

Yellow 

1368064  a  at 

Ddc 

dopa  decarboxylase  (aromatic  L-amino  acid 
decarboxylase) 

Yellow 

1368360_at 

Pig 

plasminogen 

Yellow 

1 36874  l_at 

C9 

complement  component  9 

Yellow 

1368806_at 

Seppl 

selenoprotein  P,  plasma,  1 

Yellow 

1369455_at 

Abcg5 

ATP-binding  cassette,  sub-family  G 
(WHITE),  member  5 

Yellow 

1369790_at 

Tat 

tyrosine  aminotransferase 

Yellow 

1369840_at 

UST4r 

integral  membrane  transport  protein  UST4r 

Yellow 

1370027_a_at 

LOC297568  - 
Mugl 

alpha- 1 -inhibitor  III  —  murinoglobulin  1 

Yellow 

137015  l_at 

Cpsl 

carbamoyl-phosphate  synthetase  1 , 
mitochondrial 

Yellow 

1370299_at 

Aldob 

aldolase  B,  fructose-bisphosphate 

Yellow 

1370547_at 

Pzp 

pregnancy-zone  protein 

Yellow 

1372805_at 

RGD13 10444 

LOC363015 

Green 

1373718_at 

Tubb2a 

tubulin,  beta  2a 

Green 

1373803_a_at 

Ghr 

growth  hormone  receptor 

Yellow 

1374139_at 

Cdr2 

cerebellar  degeneration-related  2 

Green 

1374650_at 

Nedd9 

neural  precursor  cell  expressed, 
developmentally  down-regulated  9 

Yellow 

1375633_at 

Clicl 

chloride  intracellular  channel  1 

Green 

1376098  a  at 

Myolg 

myosin  IG 

Green 

1377222_at 

— 

— 

Yellow 

1377461_at 

LOC680485 

hypothetical  protein  LOC680485 

Green 

1379324_at 

— 

— 

Yellow 

1379376_at 

— 

— 

Green 

1382358_at 

— 

— 

Yellow 

1383046_at 

Cfh 

complement  factor  H 

Yellow 

1383164_at 

— 

— 

Yellow 

1383169  at 

_ 

_ 

Yellow 
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APPENDIX  8  (CONTINUED):  PROBE  SETS  DIFFERENTIALLY  EXPRESSED 

BY  CHROMIUM  AND  CADMIUM 


Affymetrix  Gene  Symbol  Gene  Title  Cluster 

ID 


1385043_at 

Inadl 

InaD-like  (Drosophila) 

Yellow 

1387508_at 

Baat 

bile  acid  Coenzyme  A:  amino  acidN- 
acyltransferase  (glycine  N- 
choloyltransferase 

Yellow 

1388229_a_at 

Mugl 

murinoglobulin  1 

Y  ellow 

1388459_at 

Coll8al 

collagen,  type  XVIII,  alpha  1 

Yellow 

1389913_at 

Lrrfipl 

leucine  rich  repeat  (in  FLU)  interacting 
protein  1 

Green 

1396053_at 

Nedd9 

neural  precursor  cell  expressed, 
developmentally  down-regulated  9 

Yellow 

1384544_at 

Pon3 

paraoxonase  3 

Yellow 

1384427_at 

Mdm2 

Mdm2  p53  binding  protein  homolog 
(mouse) 

Green 
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